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ABSTRACT 


An  attack  on,  or  chemical  spill  near,  Iraq’s  oil  terminals  eould  have  disastrous 
effects  on  Iraq’s  economy.  The  impacts  from  a  drifting  mine  or  ehemical  spill  are  highly 
dependent  upon  environmental  eonditions  that  can  either  adversely  effeet  eontinued 
operations  or  hinder  the  safety  of  personnel.  Operational  Planners’  ability  to  ereate 
legitimate  seenarios  to  train  and  eombat  these  situations  is  key  to  eontinued  safe 
operation  of  the  terminals. 

To  aeeomplish  this,  planners  must  understand  the  impaets  of  using  elimatology 
versus  near  real-time  data  in  the  evaluation  of  the  seenarios.  The  near  real-time 
environmental  data  were  provided  by  the  Navy’s  oeean-atmospherie  operational  models. 
This  study  examines  the  benefits  of  knowing  the  environment  to  understand  their  impaets 
to  operations  in  the  northern  Persian  Gulf  Three  different  seenarios  were  examined  to 
understand  the  impaet  to  Oil  Terminal  operations  in  the  event  of  drifting  mines  or  a 
chemieal  spill.  The  ehemieal  spill  was  examined  from  the  outlook  of  a  major  eollision 
with  a  barge  that  spilled  either  Liquefied  Ammonia  or  Mustard  Gas.  The  Ammonia 
seenario  was  further  analyzed  in  the  atmosphere.  This  study  demonstrates  the  impaet  that 
near  real-time  environmental  eonditions  ean  have  over  elimatologieal  data  for 
Operational  Planners. 
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I. 


INTRODUCTION:  THE  PERSIAN  GULF 


A,  OVERVIEW 

The  Persian  Gulf,  also  known  as  the  Arabian  Gulf,  is  a  semienelosed  marginal  sea 
eonneeted  with  the  Indian  Ocean  through  the  Strait  of  Hormuz  with  the  major  axis 
tending  in  the  NW-SE  direction  (Figure  1).  The  Persian  Gulf  is  approximately  990  km 
long;  the  maximum  width  is  about  338  km  (Chao  et  ah,  1992).  The  estimated  surface  area 
and  volume  of  the  gulf  are  around  239,000  km^  and  8,630  km\  respectively,  which 
correspond  to  a  mean  depth  of  close  to  36  m  (Chao  et  ah,  1992).  The  maximum  depth  is 
around  100m  near  the  Strait  of  Hormuz,  with  the  Gulf  of  Oman  being  much  deeper 
(Figure  2).  The  major  axis  of  the  basin  separates  a  relatively  deeper  channel  near  Iranian 
coast  from  the  shallow  Persian  shelf  that  slops  gently  towards  the  axis  (Chao  et  al.,  1992). 


Figure  1.  General  Topography  of  region  surrounding  the  Persian  (Arabian)  Gulf  (in 

meters). 


1 


Figure  2.  General  Bathymetry  of  the  Persian  (Arabian)  Gulf  (depth  is  in  meters). 


Orographic  influences  on  the  atmospheric  circulation  are  very  significant  due  to 
the  high  mountains  bordering  much  of  the  coastline  in  this  region  (Figure  2).  The  Zagros 
Mountains  of  Iran  and  the  Jebel  al  Akhdar  range  in  northern  Oman  play  an  important  role 
in  focusing  the  winds  over  the  region.  The  Persian  Gulf  is  characterized  by  northwesterly 
or  westerly  winds  throughout  the  year  (U.S.  Hydrography  Office,  1960).  The  summer 
winds  are  mild  and  continuous.  In  contrast  the  winds  in  winter  are  often  associated  with 
synoptic  weather  systems,  and  breaks  out  suddenly  and  violently  (Reynolds,  1993). 

Due  to  its  shallow  nature,  the  Persian  Gulf  appears  to  be  influenced  by  winds  and 
surface  thermohaline  fluxes.  The  persistent  southward  wind  stress,  at  least  in  the  northern 
half  of  the  Gulf,  appears  to  set  up  coastal  current  regimes  along  both  the  Saudi 
(downwelling)  and  Iranian  (upwelling)  coasts.  A  persistent  thermal  front  across  the 
Persian  Gulf  about  the  latitude  of  Qatar  appears  related  to  the  thermohaline  exchange 
through  the  Hormuz  Strait  and  splits  the  Persian  Gulf  into  two  regimes. 
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The  eastern  or  southern  regime  is  a  dynamic  regime.  It  lies  on  the  border  between 
tropical  and  extratropical  regimes  and  feels  the  influences  from  both.  The  winds  and 
currents  in  the  Gulf  of  Oman  are  influenced  by  the  southwest  monsoon.  This  influence 
continues  through  the  Straights  of  Hormuz  and  into  the  eastern  regime.  The  northern,  or 
western,  regime  is  characterized  by  low  energy  throughout.  This  includes  both  currents 
and  winds.  The  winds  are  predominantly  extratropical  in  nature,  lying  on  the  fringe  or 
extratropical  influences.  The  currents  tend  to  be  extremely  low.  Influences  from  the 
eastern  regime  occasionally  reach  the  northern  regime  which  increases  the  energy  of  the 
northern  regime.  Otherwise  the  flow  remains  low  energy. 

B,  ECONOMIC  AND  POLITICAL  CONCERNS 

The  Persian  Gulf  region  continues  to  be  of  paramount  importance  to  the  world 
community  for  a  variety  of  economic  and  political  reasons.  The  immense  petroleum 
resources  and  their  maritime  export  from  the  region  and  the  complex  evolution  and 
conflict  of  competing  political  entities  are  constantly  demanding  attention. 

There  have  been,  and  continue  to  be,  significant  territorial  disputes  between 
Persian  Gulf  countries.  Two  of  the  major  disputes  in  recent  history  involve  the  Iraqi 
invasion  of  Kuwait  in  August  1990  and  the  Iran-Iraq  War  from  1980  to  1988.  Unlike 
Europe,  the  boundaries  of  the  countries  were  not  derived  from  language  or  cultural 
influences.  In  some  cases  the  boundaries  were  decided  by  countries  outside  the  region 
after  World  War  II.  This  compounds  the  problem,  creating  numerous  religious  and 
ethnical  backgrounds  within  each  country  that  do  not  necessarily  want  to  coexist. 

According  to  the  Energy  Information  Administration  International  Energy 
Outlook  2003  (EIA,  2003)  the  Persian  Gulf  countries  produced  27%  of  the  world  oil, 
while  holding  57%  of  the  world’s  proven  crude  oil  reserves.  Besides  oil,  the  Persian  Gulf 
region  also  has  huge  reserves  of  natural  gas,  accounting  for  45%  of  the  total  proven  world 
gas  reserves.  The  Energy  Information  Administration  International  Energy  Outlook  2003 
predicts  the  Persian  Gulf  oil  production  capacity  would  increase  to  33%  of  the  world  total 
by  2020. 
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With  close  to  a  third  of  the  world  petroleum  eoming  from  the  Persian  Gulf 
countries,  the  vast  majority  of  it  is  exported  via  tanker  through  the  Strait  of  Hormuz. 
According  to  Energy  Information  Administration  International  Energy  Outlook  2003 
(EIA,  2003),  the  amount  transiting  through  this  strait  was  close  to  90%  of  all  oil  coming 
from  the  region.  Of  primary  importance  to  this  study  is  Iraq’s  contribution. 

Iraq's  economy  is  dominated  by  the  oil  seetor,  whieh  has  traditionally  provided 
about  95%  of  foreign  exchange  earnings.  In  the  Persian  Gulf,  Iraq  has  three  tanker 
terminals:  Bakr  A1  Basrah  (ABOT),  Khor  al-Amaya  (KAAOT),  and  Khor  az-Zubair 
(which  mainly  handle  dry  goods  and  minimal  oil  volumes).  ABOT,  formerly  known  as 
Mina  al,  is  Iraq's  largest  oil  terminal,  with  two  pipelines  (48-inch  and  41 -inch),  plus  four 
400,000-bbl/d  capacity  berths  capable  of  handling  very  large  erude  carriers  (EIA,  2006). 
ABOT  is  loeated  on  the  northern  tip  of  the  Persian  Gulf,  at 
approximately 29. 75°  A  48.83°E’,  with  KAAOT  loeated  just  to  the  northwest,  at 
approximately  29.19° N  48.80°E'  (BP  Port  Eist,  2007).  The  loeations  of  the  oil  terminals 
are  plotted  in  Eigure  3. 


Eigure  3.  Plot  of  Iraqi  Oil  Terminals  ABOT,  five  -pointed  white  star,  and  KAAOT, 
six-pointed  yellow  star  (depth  is  in  meters). 
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C.  SCENARIO  SELECTION  CONCEPT 

The  execution  of  any  naval  operation  can  be  hindered  by  numerous  environmental 
factors.  Proper  planning  for  these  operations  includes  forecasting  the  environment  of  the 
area  of  operations.  Where  drifting  mines  or  chemical  spills  are  concerned  the 
environment  plays  an  important  role  in  determining  spill  extent  or  mine 
avoidance/placement. 

As  the  warfare  environment  continues  to  shift  from  the  open  ocean  to  littoral, 
environmental  prediction  and  modeling  efforts  of  the  shallow  water  surroundings  need  to 
be  understood  and  utilized  in  order  to  support  these  operations.  The  difference  between 
climatological  data  and  observations  can  be  dramatic.  Knowledge  of  the  environment  and 
the  ability  to  predict  it  can  determine  the  success  of  an  operation  in  many  instances.  The 
scenarios  considered  in  this  thesis  are  chosen  to  emphasize  the  importance  of  knowing, 
and  predicting,  environmental  factors  to  maximize  the  operation’s  success  rate. 

High  resolution  full  physical  models  are  needed  to  forecast  the  surface  winds  and 
currents.  After  the  high-resolution  winds  and  currents  are  obtained,  the  transport  models 
and  chemical  models  are  utilized  to  forecast  the  mine  drift  and  chemical  spill.  Two 
physical  and  two  chemical  models  were  utilized  for  this  thesis.  The  two  physical  models, 
(1)  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  (COAMPS),  and  (2)  the 
Shallow  Water  Assimilation  Forecast  System  (SWAPS),  are  used  to  get  the  physical 
environment  such  as  the  surface  winds  and  surface  currents,  respectively.  The  outputs  of 
the  physical  models  are  then  used  as  inputs  into  the  two  chemical  models:  (1)  the  Hazard 
Prediction  and  Assessment  Capability  (HP AC)  for  the  atmosphere  and  (2)  CHEMMAP^'^ 
for  the  ocean.  The  outputs  were  then  analyzed  and  applied  to  one  of  several  possible 
naval  operations. 

Statistics  and  complex  Empirical  Orthogonal  Eunction  (EOE)  analyses  are  used  to 
classify  the  physical  environment  from  the  COAMPS  and  SWAPS  such  as  the  winds  and 
currents  in  the  Persian  Gulf  during  the  period  from  the  beginning  of  Pebruary  through  the 
end  of  July,  2006.  The  purpose  is  to  examine  the  extreme  environmental  effects  on 
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different  possible  seenarios  in  the  northern  seetion  of  the  Persian  Gulf.  The  seenarios 
inelude  a  drift  seenario  and  different  ehemieal  spills,  either  aeeidental  or  planned. 

The  drifting  mine  seenario  looks  at  the  possibility  of  floating  mines  affeeting  oil 
terminal  operations.  This  seenario  ineludes  a  generie  one  meter  diameter  mine  with  a 
density  of  Xfygjcm  .  This  equates  to  eaeh  mine  weighing  approximately  one  ton.  This 
generie  mine  will  then  be  set  adrift  in  various  spatial  and  temporal  loeations  throughout 
the  northern  Persian  Gulf  The  seenario  will  also  examine  the  ability  to  utilize 
elimatology  as  a  valid  first  guess. 

The  other  portion  of  the  study  ineludes  two  different  ehemieal  spills  loeated  near 
the  oil  terminals.  These  seenarios  will  be  run  during  times  that  meet  eertain  eriteria.  This 
will  in  turn  ensure  worst  ease  seenarios  are  examined.  The  eriteria  will  inelude;  (1)  low 
winds  and  weak  eurrents,  (2)  low  winds  and  strong  eurrents,  (3)  high  winds  and  weak 
eurrents,  (4)  strong  winds  and  strong  eurrents,  and  (5)  speeifie  events  sueh  as  frontal 
passage.  These  events  will  be  eompared  to  elimatologieal  data  to  determine  how  mueh 
value  ean  be  added  over  elimatology. 

All  of  these  seenarios  will  be  examined  through  the  eoneept  of  what  and  how 
mueh  value  ean  be  added  by  near  real-time  data  versus  elimatologieal  or  no  data.  This 
eoneept  is  an  underlying  eoneept  behind  ‘Battlespaee  on  Demand’.  The  purpose  of 
‘Battlespaee  on  Demand’  is  to  effeetively  supply  eombatant  eommanders  with 
information  superiority  over  the  entire  area  of  operation.  This  eoneept  is  then  broken 
down  into  three  tiers,  whieh  inelude:  (1)  Environment  Layer,  (2)  Performanee  Layer,  and 
(3)  Deeision  Layer.  Ligure  4  illustrates  the  Three  Tiers  of  Battlespaee  on  Demand.  This 
thesis  utilized  this  eoneept  in  the  analysis  of  the  data  with  respeet  to  the  seenarios.  The 
ability  to  add  value  by  using  real-time  or  elimatologieal  data  was  of  interest  when 
analyzing  the  seenarios. 
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Battlespace  on  Demand 

The  Three  Tiers 


Tier  3  -  the  Decision  Layer 

•  Options  /  Courses  of  Action 

•  Search  Patterns 

•  Asset  Allocation  /  Timing 

•  Quantify  Risk 


Tier  2  -  the 
Performance  Layer 


Tier  1  -  the 
(forecast) 
Enviromneni 
Layer 


Initial  and  Boundary  Conditions 


- United  States  Fleet  Forces 


Operational  Readiness.  EffecWeness.  Primacy  — 


Figure  4.  Three  Tiers  of  Battlespace  on  Demand  (From  Capt.  Dave  Titley,  USN, 

2007). 
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II.  ENVIRONMENTAL  FACTORS  FROM  PREVIOUS  STUDIES 


A.  CURRENTS 

The  basic  features  of  the  Persian  Gulf  can  be  divided  into  a  northern  and  southern 
or  eastern  regime  (Figure  5).  The  primary  focus  of  this  research  is  the  northern  regime. 
The  northern  regime  is  dominated  by  wind  forcing  to  the  south  along  the  axis  of  the  Gulf 
and  the  riverine  input  at  the  Gulfs  head.  The  wind-driven  response  of  the  Gulf  appears  to 
be  the  typical  adjustment  of  the  pressure  field  such  as  to  produce  a  down-wind  flow,  i.e., 
there  is  downwelling  on  the  western  coast  and  upwelling  on  the  coast  of  Iran,  and 
evidence  for  a  southeastward  flowing  coastal  current  along  both  the  northern  and 
southern  coasts  (Reynolds,  1993). 


Figure  5.  Sketch  of  surface  currents  and  circulation  processes  in  the  Persian  Gulf  and 
western  Gulf  of  Oman  (From  Reynolds,  1993). 

The  flow  along  the  Kuwait  and  Saudi  coast  is  augmented  by  the  freshwater  input 
from  the  north  which  forms  a  riverine  plume.  The  river  inflows  are  approximately  split 
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between  the  flow  out  of  the  Shalat  Ariabi  (Tigris  and  Euphrates)  and  rivers  flowing  out  of 
the  highland  of  Iran  (the  Hendijan,  Hilleh,  and  Mand).  Today  the  flow  of  the  Shalat 
Ariabi  is  mueh  smaller  than  it  onee  was  beeause  of  massive  dam  projeets  in  Turkey.  It  is 
not  yet  elear  what  ehanges  this  deeline  in  freshwater  input  has  made.  The  eenter  of  the 
northern  Gulf  appears  to  be  fairly  stagnant  (Reynolds,  1993).  The  southern  end  of  this 
regime  eorresponds  roughly  to  the  longitude  of  Qatar  and  Bahrain,  although  the 
termination  of  the  northern  eireulation  is  poorly  understood.  The  flow  along  the  Iran 
eoast  seems  to  eontinue  into  the  southeastern  basin  as  a  tightly  trapped  eoastal  eurrent 
extending  perhaps  as  far  as  the  Strait  of  Hormuz. 

The  northern  Persian  Gulf  is  separated  from  the  southern  regime  by  a  front  off 
Qatar  (Figure  5),  whieh  is  most  intense  in  summer  and  weakest,  at  least  in  sea  surfaee 
temperature,  in  the  late  winter  and  spring  (Reynolds,  1993).  The  loeation  of  this  front 
appears  in  both  elimatologieal  hydrographie  data  and  remotely  sensed  SST  to  be  tied  to 
the  penetration  of  fresh  inflow  into  the  Gulf  from  the  Strait  of  Hormuz. 

B,  WINDS 

The  Gulf  is  loeated  between  latitudes  24°  A  and 30°  A  ,  in  whieh  are  loeated  most 
of  the  Earth's  deserts.  This  region  marks  the  boundary  between  tropieal  eireulations 
(Hadley  Cell)  and  the  synoptie  weather  systems  of  mid-latitudes.  Deseending  dry  air  in 
these  latitudes  produees  elear  skies  and  arid  eonditions. 

In  the  north,  the  loeal  elimate  is  influeneed  by  orography.  The  Taurus  and  Pontie 
mountains  of  Turkey,  and  the  Zagros  Mountains  of  Iran  together  with  the  Tigris- 
Euphrates  Valley,  form  a  northwest-southeast  axis  that  strongly  influenees  the  traeks  of 
extra-tropieal  storms  to  a  southeast  direetion.  As  part  of  an  analysis  from  1978  to  1998  of 
mean  synoptie  eonditions,  it  was  founds  that  winds  in  the  northern  region  predominantly 
follow  the  major  axis  of  the  Persian  Gulf.  This  is  elearly  illustrated  in  Kuwait 
International  Airport  (K)  portion  of  Figure  6  (Zhu  and  Atkinson,  2004).  Figure  7 
illustrates  the  same  eoneept;  however  the  mean  and  peak  speeds  are  also  shown.  Form 
Figure  7,  the  elimatologieal  mean  ean  be  seen  to  vary  from  10  to  15  knots,  or  around  5  to 
8  m/s. 
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The  most  well-known,  and  notorious,  weather  phenomenon  in  the  Gulf  is  the 
Shamal,  a  NW  wind  whieh  oeeurs  year  round  (Perrone,  1979).  Shamal  is  an  Arabie  word 
meaning  north.  The  Winter  Shamal  is  a  wind  that  sets  in  with  great  abruptness  and  foree, 
and  is  related  to  synoptie  weather  systems  to  the  NW.  It  seldom  exeeeds  10  m/s  (<5% 
frequeney)  but  lasts  several  days  (Reynolds,  1993).  The  summer  Shamal  is  praetieally 
eontinuous  from  early  June  through  July.  It  is  assoeiated  with  the  relative  strengths  of  the 
Indian  and  Arabian  Peninsula  thermal  lows. 


E 


E 


Figure  6.  Frequeney  of  the  hourly  wind  direetion  and  ealms  (in  eenter)  observed  at 
Kuwait  International  Airport  (K)  in  2002:  (a)  January;  (b)  April;  (c)  July;  (d)  Oetober 

(After  Zhu  and  Atkinson,  2004). 
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KU\A)AIT  INTL  (CIV/WIL)  KU  AFCCC  (ARAIA) 

Block  Station:  +06820  ICAO:  OKBK  Lot:  29,22  N  Lon:  47.98  E  151  Patton  Ave,  Rm  120 

Dota  Derived  from  AFCCC  Operational  Climatic  Dato  Summary  (OCDS)  Asheville,  NC  28801-5002 


Figure  7.  Climatological  summary  of  winds  at  Kuwait  International  airport.  Scale  on 
the  right  is  for  peak  winds  and  range  between  30  and  70  knots  (From  AFCCC 
https://notus2.afccc.afmil/SCISPublic/,  2007). 


The  winter  Shamal  brings  some  of  the  strongest  winds  and  highest  seas  of  the 
season  to  the  Gulf  region  (Reynolds,  1993).  Winds  in  the  area  ahead  of  an  approaching 
cold  front  blow  from  the  southeast.  These  winds,  called  Kaus  in  Arabic  or  Shakki  in 
Arabian,  slowly  increase  in  intensity  as  the  front  approaches.  Winds  may  reach  gale  force 
before  the  passage  of  the  front  and  the  onset  of  the  Shamal.  Due  to  the  channeling  of  the 
low-level  air  flow  by  the  Zagros  Mountains  of  western  Iran,  the  strongest  of  the  southerly 
winds  occur  on  the  eastern  seaboard.  The  Shamal  usually  occurs  first  in  the  northwest 
and  then  spreads  south. 

Zhu  and  Atkinson  (2004)  completed  a  study  of  observed  and  modeled 
climatology  in  the  Persian  Gulf.  Part  of  that  study  included  Kuwait  International  Airport 
(KWI),  located  at29.23°  A  48.97°£' .  In  January  (Figure  6a),  wind  direction  at  KWI  was 
predominantly  northwesterly.  In  April  (Figure  6b),  at  KWI,  southeasterly  and  easterly 
directions  were  now  of  significant  frequency.  In  July  (Figure  6c),  northwesterly  and 
northerly  directions  again  dominated  the  area  near  KWI.  By  October  (Figure  6d),  easterly 
and  southeasterly  directions  were  more  evident  at  KWI.  Northwesterly  flow  was 
dominant  for  most  of  the  year  over  the  entire  Gulf,  except  the  eastern  area  near  Bandar 
Abbas.  Over  the  west  coast  of  the  Gulf  the  picture  was  not  clear  cut.  At  KWI  in  January 
the  direction  of  the  mean  wind  was  clearly  from  the  northwest,  but  on  some  occasions  the 
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wind  blew  from  the  east  between  1400  LT  (1100  UTC)  and  2000  LT  (1700  UTC).  In 
April,  the  direetion  of  the  daytime  mean  wind  was  from  the  east  and  southeast,  i.e.  from 
the  sea.  In  July,  the  mean  direetion  was  again  from  the  northwest  throughout  the  day.  The 
Oetober  pattern  was  rather  similar  to  that  in  April. 

A  strong  sea  breeze  oeeurs  along  the  entire  eoastline,  espeeially  along  the  Arabian 
Peninsula.  Driven  by  the  intense  temperature  differenee  between  the  land  and  water 
surfaees,  the  sea  breeze  eireulation  adds  a  landward  eomponent  to  all  winds.  The  effeet  of 
these  winds  is  to  drive  surfaee  pollutants  (oil  and  others)  to  the  beaeh  mueh  faster  than 
they  would  move  otherwise.  Examples  of  on-  or  off-shore  winds  as  high  as  15  m/s  are 
evident  in  eoastal  meteorologieal  stations  (Reynolds,  1993). 

C.  TIDES 

The  tides  in  the  Gulf  eo-oseillate  with  those  in  the  narrow  Strait  of  Hormuz, 
whieh  opens  into  the  deep  Gulf  of  Oman  (Figure  8).  The  tides  in  the  Gulf  are  eomplex 
standing  waves  and  the  dominant  pattern  varies  from  being  primarily  semi-diurnal  to 
diurnal.  The  tidal  range  is  large,  with  values  greater  than  1 .0  m  everywhere  (Reynolds, 
1993).  The  dimensions  of  the  Gulf  are  sueh  that  resonanee  amplifieation  of  the  tides  ean 
oeeur  and  the  result  is  that  the  semi-diurnal  eonstituents  have  two  amphidromie  points,  in 
the  northwest  and  southeast  ends,  and  the  diurnal  eonstituents  have  a  single  amphidromie 
point  in  the  eenter  near  Bahrain. 

Tides  are  important  for  stirring  and  mixing  waters  vertieally  and  on  a  horizontal 
seale  of  10  km,  but  tides  are  not  an  important  eontribution  to  the  residual  eireulation  of 
the  Gulf  (Reynolds,  1993).  Tidal  eurrent  (averaging  longer  than  a  day)  has  negligible 
residual  energy,  and  as  a  result,  basin-seale  adveetion  from  tides  is  not  eonsidered  by 
most  oil  trajeetory  models  or  general  eireulation  models.  Tides  are  important  on  smaller 
seales  of  horizontal  length  (<  10  km)  and  time  (<24  hr). 
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Figure  8.  The  M^,S2,  and  tidal  eonstituents  in  the  Gulf  (After:  Reynolds.,  1993). 

Tide  heights  are  shown  in  meters. 


D,  CURRENTS:  TIDAL,  WIND-DRIVEN,  AND  DENSITY  DRIVEN 

The  energy  in  the  water  motion  ean  be  related  to  three  foreing  proeesses:  tidal 
forees,  wind  forees,  and  density  differenees.  The  kinetie  energy  of  the  water  veloeity  ean 
be  partitioned  among  the  three  terms  approximately  as  100,  10,  and  1  respeetively 
(Reynolds,  1993).  Eaeh  of  the  different  eurrents  has  a  different  sealing  time:  tides  vary 
over  a  few  hours  at  diurnal  or  semi-diurnal  periods,  wind-driven  eurrents  develop  and 
subside  over  a  few  days;  and  density-driven  eurrents  take  weeks  to  ehange  in  response  to 
seasonal  foreing.  These  three  eurrent  eomponents  ean  be  extraeted  from  measurements, 
but  the  seleetion  of  averaging  times  and  filter  types  affeets  the  eomputed  proportions 
eonsiderably  and  must  by  used  with  eare.  Beeause  of  their  eonsistent  direetion,  density 
eurrents  are  important  in  the  distribution  and  removal  of  pollutants  from  the  Gulf  region. 
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The  main  features  of  the  residual  eireulation  in  the  Gulf  are  (Hunter,  1983):  (1) 
high-  and  low-salinity  water  exehange  in  the  Strait  of  Hormuz;  (2)  density-dominated 
eireulation  in  the  eentral  and  southern  regions;  (3)  frietional-balaneed,  wind-dominated 
eireulation  in  the  NW  region;  and  (4)  evaporation-indueed  bottom  flow.  The  aetual 
pattern  of  eireulation  is  complex. 

An  important,  secondary  current  in  the  Gulf  is  a  coastal,  reverse  circulation  along 
the  Iranian  coast.  Driven  by  the  density  differences  that  result  from  river  runoff,  the 
coastal  currents  along  the  Iranian  coast  are  southerly,  against  the  inflow  water.  The 
reverse  currents  are  evident  in  shuttle  photographs  and  Advance  Very  High  Resolution 
Radiometer  (AVHRR)  images  (Reynolds,  1993).  The  interface  between  the  two  currents 
is  a  region  of  horizontal  current  shear  which  often  shows  evidence  of  wavelike  shear 
instability. 

E,  MIXING  PROCESSES  AND  VERTICAL  STRUCTURE, 

There  are  three,  primary  contributors  to  mixing  to  the  water  column:  (1)  tides,  (2) 
winds  and  waves,  and  (3)  evaporation.  Internal  waves  at  density  interfaces,  and 
topographic  features  such  as  islands  also  contribute  to  mixing,  but  their  contribution  is 
unknown  and  probably  small  (Reynolds,  1993).  The  strong  tidal  currents  create  a  bottom, 
mechanical,  turbulent  friction  layer.  The  enhanced  turbulent  mixing  homogenizes  the 
water  column  from  the  bottom  up,  and  a  bottom  mixed  layer  is  evident  in  all 
hydrographic  sections.  The  other  two  mechanical  mixing  processes  are  strongest  at  the 
surface.  Wind  friction  creates  a  surface  shear  layer  and  also  creates  surface  wave  fields, 
both  of  which  form  a  surface  mixed  layer.  Evaporation  of  fresh  water  enhances  mixing 
by  increasing  the  salinity  and  density  of  the  surface  water,  thus  reducing  surface  stability. 
In  the  winter,  combined  cooling  and  evaporation  increase  the  density  sufficiently  to 
overturn  the  water  column  and,  north  of  Qatar,  create  a  well-mixed  water  column. 

Mixing  processes  are  opposed  by  the  stabilizing  effects  of  surface  warming.  As 
the  temperature  increases,  the  water  density  decreases.  Countering  this  stabilization, 
evaporation  increases  surface  salinity,  which  increases  density  and  creates  an  unstable 
overturn  to  some  lower  depth.  The  resulting  surface  mixed  layer  and  lower  thermocline 
reach  about  20  m  depth  (Reynolds,  1993). 
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III.  OCEAN-ATMOSPHERIC  PHYSICAL-CHEMICAL  MODELS 


This  thesis  utilizes  two  physieal  models  and  two  ehemieal  models.  Among  them, 
COAMPS  and  HP  AC  are  the  physieal  and  ehemieal  models  for  the  atmosphere;  and 
SWAPS  and  CHEMMAP™  (developed  by  the  Applied  Seienee  Assoeiates  (ASA),  Ine) 
are  the  physieal  and  ehemieal  models  for  the  oeean.  The  model  output  of  COAMPS  sueh 
as  the  winds  is  the  input  for  the  atmospherie  ehemieal  model  (HP AC)  and  the  oeeanie 
physieal  and  ehemieal  models:  SWAPS  and  CHEMMAP.  To  solely  study  the  effeets  of 
atmospherie  dispersion  HP  AC  was  utilized.  To  understand  the  impaets  of  surfaee  eurrents 
and  winds  on  mine  drifting  and  ehemieal  spill,  the  Eagrangian  drift  model  and 
CHEMMAP  were  employed.  The  following  seetion  will  deseribe  the  model 
eharaeteristies  and  uses. 

A,  COAMPS 

The  following  overview  of  COAMPS  is  taken  direetly  from  Chen  et  al.  (2003). 
The  atmospherie  eomponent  of  COAMPS,  the  portion  used  for  this  analysis,  ean  be  used 
for  real-data  applieations.  Por  these  applieations,  the  COAMPS  analysis  ean  use  global 
fields  from  the  Navy  Operational  Global  Atmospherie  Predietion  System  (NOGAPS) 
foreeast  as  the  first-guess.  Observations  from  aireraft,  rawinsondes,  ships,  and  satellites 
are  blended  with  the  first-guess  fields  to  generate  the  eurrent  analysis.  The  atmospherie 
model  uses  nested  grids  to  aehieve  high  resolution  for  a  given  area;  it  eontains 
parameterizations  for  subgrid  seale  mixing,  eumulus  parameterization,  radiation,  and 
explieit  moist  physios. 

The  atmospherie  portion  of  the  COAMPS  models  is  oomprised  of  the 
nonhydrostatio,  fully  oompressible  equations  of  motion  following  Klemp  and 
Wilhelmson  (1978).  The  adiabatio  equations  are  developed  using  the  equation  of  state: 

p^pRJ^,  (3.1) 

where  p  is  the  pressure,  p  is  the  density,  is  the  dry  gas  oonstant,  the  virtual 
temperature,  L  is 
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t;  =  r  (1.0+ 0.608^  J, 


(3.2) 


where  T  is  the  temperature  and  is  the  speeifie  humidity,  and  the  Exner  funetion, 


(3.3) 


where  ep  is  the  speeifie  heat  at  eonstant  pressure,  and  p^g  the  referenee  pressure. 

The  transformation  of  the  vertieal  eoordinate  following  Gal-Chen  and  Somerville 
(1975)  is  applied  to  map  the  lowest  eoordinate  surfaee  to  an  irregular  lower  boundary 


where  is  the  depth  of  the  model  domain  and  is  the  height  of  the  topography. 

Both  the  horizontal  and  vertieal  grids  in  the  foreeast  model  are  staggered.  The  horizontal 
grid  uses  the  Arakawa  and  Lamb  (1977)  seheme  C  staggering  depieted  in  the  Figure  9. 
The  sealar  variables  (denoted  ;r)  are  defined  in  the  eenter  of  the  model  grid  box: 
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Figure  9.  Diagram  of  Arakawa  scheme  C  (From:  Comet  Module 
http://www.comet.ucar.edu/,  2007). 


The  computational  grid  features:  (1)  the  w-component  are  one-half  grid  intervals 
between  the  mass  points  in  the  x  direction,  (2)  the  v-component  are  one-half  grid  intervals 
between  the  mass  points  in  the  y  direction,  and  (3)  the  w-component,  are  coincident  with 
the  mass  variables.  In  general,  all  derivatives  are  computed  to  second-order  accuracy. 
Exceptions  are  the  horizontal  diffusion  and  an  option  for  fourth-order  accurate  horizontal 
advection.  By  using  fourth  order  accurate  horizontal  diffusion,  the  damping  is  much  more 
specific  to  the  removal  of  high-frequency  modes. 

Solving  the  fully  compressible,  nonhydrostatic  equations  explicitly  is  extremely 
computational  expensive  because  of  the  presence  of  sound  waves  that  severely  limit  the 
time  step  required  to  maintain  computational  stability.  One  approach  for  addressing  this 
problem  is  to  treat  the  sound  wave  modes  separately  on  a  small  time  step  following 
Klemp  and  Wilhelmson  (1978)  and  Skamarock  and  Klemp  (1992). 

The  equations  are  solved  using  the  centered-in-time  or  leapfrog  scheme  (Haltiner 
and  Williams  1980).  Typically,  the  time  splitting  due  to  the  leapfrog  scheme  is  not 
problematic  (e.g.,  Klemp  and  Wilhelmson  1978).  However,  a  Robert  (1966)  time 
smoother  is  applied  to  assure  that  any  tendencies  that  may  tend  to  decouple  the  odd  and 
even  time  steps  are  stable.  For  any  variable  (j) ,  the  Robert  time  filter  is  applied  as 

<!)’  =f  -If 

The  first  equation  corresponds  to  the  leapfrog  step  for 


with  the  asterisk  corresponding  to  provisional  terms  that  have  not  yet  been  smoothed 
through  application  of  the  second  step.  The  net  effect  is  to  produce  strong  damping  of  the 
computational  mode,  while  the  physical  mode  is  generally  not  affected  (Asselin  1972). 

For  typical  COAMPS  applications,  a  is  0.2. 
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COAMPS  uses  a  level  2.5  seheme  (Mellor  and  Yamada  1982)  that  solves  both  a 
prognostie  equation  for  turbulent  kinetie  energy  (IKE)  and  diagnostie  equations  for 
seeond-moment  quantities  sueh  as  primarily  fluxes  of  heat,  moisture,  and  momentum. 
Sehemes  that  feature  a  prognostie  equation  for  TKE  and  diagnostie  equations  for  other 
quantities  are  referred  to  as  “1.5  order  elosure.”  Thus,  COAMPS  eontains  a  1.5  order 
elosure,  level  2.5  seheme. 

The  1.5  order  elosure,  level  2.5  seheme  eomputes  all  of  the  fields  neeessary  for 
solving  the  TKE  equation,  ineluding  boundary  layer  depth,  turbulent  mixing  length,  flux, 
Riehardson  number,  and  eddy  eoeffieients.  In  addition,  the  seheme  ineorporates  the 
influenee  of  boundary  layer  eloudiness. 

The  surfaee  layer  parameterization  follows  the  Eouis  (1979)  seheme,  whieh  uses 
polynomial  funetions  of  the  bulk  Riehardson  number  to  direetly  eompute  surfaee  sensible 
heat  flux,  surfaee  latent  heat  flux,  and  surfaee  drag.  The  bulk  Riehardson  number  is 
defined  as 


Rig 


(3.6) 


where  g  is  the  aeeeleration  due  to  gravity,  z  is  the  referenee  elevation  (equal  to  10  m  in 
COAMPS),  is  the  air-sea  temperature  differenee,  u  is  the  wind  speed  at  the  referenee 
elevation,  and  0  is  the  mean  potential  temperature  over  the  depth  of  the  surfaee  layer. 

Surfaee  roughness  (Eairall  et  al.  1996)  is  obtained  by 

^0=^0— +  C  —  ’  (3-7) 

g 


where  e^,  is  the  Chamoek  eonstant,  u*  is  the  frietion  veloeity,  g  is  the  aeeeleration  due  to 
gravity,  e^,  is  a  eonstant,  and  v  is  the  moleeular  viseosity.  The  first  term  is  the  Charnoek 
relation  that  aeeounts  for  high  wind  (aerodynamieally  rough  eonditions),  while  the 
seeond  term  aeeounts  for  low  wind  (aerodynamieally  smooth  eonditions). 
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Although  the  equations  are  solved  on  a  staggered,  seheme  C  grid,  the  COAMPS 
analysis  is  performed  on  the  Arakawa-Lamb  seheme  A  grid  (i.e.,  no  grid  staggering).  The 
bieubie  spline  interpolation  is  used  to  interpolate  the  analyzed  fields  to  the  C  grid  within 
the  foreeast  model  eode. 

The  COAMPS  analysis  is  based  on  the  multivariate  optimum  interpolation 
(MVOI)  analysis  seheme  deseribed  in  Goerss  and  Phoebus  (1992)  and  Barker  (1992). 
The  MVOI  teehnique  uses  observational  data  to  eompute  inerements  for  the  first-guess 
fields.  If  eomputed  over  a  statistieally  signifieant  number  of  oases,  it  minimizes  the  mean 
squared  error  of  the  analysis.  The  analysis  variables  for  the  MVOI  are  geopotential 
height,  and  the  u  and  v  wind  oomponents.  Finally,  the  first-guess  fields  are  adjusted  based 
on  observational  data  via  a  MVOI  analysis.  Observational  data  inolude  the  following  data 
types:  (1)  Radiosonde,  (2)  Pibal,  (3)  Surfaoe  land,  (4)  Surfaoe  marine,  (5)  Airoraft,  (6) 
Satellites  (inoluding  SSM/I,  Soatterometer,  Sea  Surfaoe  Temperature,  and  QUIKSoat), 
and  (7)  Synthetio  observations  from  NOGAPS.  The  outoff  for  data  used  is  ±3: 00  hours 
from  the  analysis  time.  Conventional  data  are  subjeoted  to  quality  oontrol  (Baker,  1992 
and  1994)  inoluding  gross  oheok  error  oheoking  and  oomplex  quality  oontrol  of 
radiosonde  observations  (Gandin,  1988).  Quality  oontrol  of  airoraft  data  inoludes 
sophistioated  flight  traok  oheoking  and  oharaoteristio  error  deteotion  (Pauley  2001). 
Within  the  MVOI  itself,  both  satellite  and  oonventional  data  are  further  oheoked  for 
quality  and  for  oonsistenoy  with  neighboring  observations  and  the  model  short-term 
foreoast,  whioh  is  used  as  a  "first  guess"  for  the  analysis. 

After  going  through  quality  oontrol,  the  time/spaoe  adjusted  model  foreoast 
variables  are  oompared  to  the  observations  to  determine  differenoes  to  be  input  into  the 
analysis  program.  The  differenoes  are  then  weighted  by  a  faotor  based  on  the  reliability  of 
the  observations  (instrument  error  oharaoteristios,  for  example)  relative  to  that  of  the  first 
guess.  A  mathematioal  prooedure  is  then  employed  to  minimize  the  analysis  error,  based 
on  the  assumed  observation  and  model  foreoast  error  oharaoteristios.  The  resulting 
analysis  inerements  are  then  applied  to  the  first  guess,  resulting  in  a  new  analysis. 
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B,  SWAFS 

SWAPS  is  a  coastal  ocean  forecast  system  in  operational  use  at  the  Naval 
Oeeanographie  Offiee  (Clifford,  et  ah,  1994).  The  modeling  system  is  built  based  on  the 
Prineeton  Oeean  Model  (POM)  with  the  main  differenees  in  data  assimilation. 

The  following  overview  of  POM  is  taken  direetly  from  Mellor  (2004).  The 
prineipal  attributes  of  the  model  are  as  follows:  (1)  It  eontains  an  imbedded  seeond 
moment  turbulenee  elosure  sub-model  to  provide  vertieal  mixing  eoefficients,  (2)  It  is  a 
sigma  eoordinate  model  in  that  the  vertieal  eoordinate  is  sealed  on  the  water  eolumn 
depth,  (3)  The  horizontal  grid  uses  eurvilinear  orthogonal  eoordinates  and  an  "Arakawa 
C"  differeneing  seheme,  (4)  The  horizontal  time  differeneing  is  explieit  whereas  the 
vertieal  differeneing  is  implieit.  The  latter  eliminates  time  eonstraints  for  the  vertieal 
eoordinate  and  permits  the  use  of  fine  vertieal  resolution  in  the  surfaee  and  bottom 
boundary  layers.  (5)  Complete  thermodynamies  have  been  implemented. 

The  turbulenee  elosure  sub-model  is  one  that  Mellor  and  then  was  signilieantly 
advaneed  in  eollaboration  with  Tetsuji  Yamada  (Mellor  and  Yamada,  1974;  Mellor  and 
Yamada,  1982).  Here,  the  Level  2.5  model  is  used  together  with  a  prognostie  equation  for 
the  turbulenee  maeroseale. 

By  and  large,  the  turbulenee  model  seems  to  do  a  fair  job  simulating  mixed  layer 
dynamies  although  there  have  been  indieations  that  ealeulated  mixed  layer  depths  are  a 
bit  too  shallow  (Martin,  1985).  Also,  wind  foreing  may  be  spatially  and  temporally 
smoothed.  It  is  known  that  the  latter  proeess  will  reduee  mixed  layer  thieknesses  (Klein, 
1980). 

The  basie  equations  have  been  east  in  a  bottom  following,  sigma  eoordinate 
system  whieh  is  illustrated  in  Figure  10.  The  reader  is  referred  to  Phillips  (1957),  and 
Blumberg  and  Mellor  (1980,  1987)  for  a  derivation  of  the  sigma  eoordinate  equations 
whieh  are  based  on  the  transformation, 

*  *  ^  /7  * 

=  — (3.8) 
H  ^-7] 
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where  x,y,z  are  the  eonventional  Cartesian  eoordinates;  where  //(x,j)is  the  bottom 
topography  and  ri{x,y,t)  is  the  surfaee  elevation.  Thus,  cr  ranges  from  cr  =  0  at  z  =  ;/  to 

cr  =  -l  ai  z  =  H .  The  sigma  eoordinate  system  is  probably  a  neeessary  attribute  in 
dealing  with  signifieant  topographieal  variability  sueh  as  that  eneountered  in  estuaries  or 
over  eontinental  shelf  breaks  and  slopes. 


Let  {U,  F)  be  the  horizontal  veloeity  eomponents;  p  be  the  density;  (Km,  Kh)  be 
the  eddy  viseosity  and  thermal  diffusivity.  After  eonversion  to  sigma  eoordinates  and 
deletion  of  the  asterisks,  the  eontinuity  equation  is  given  by 


dDU  SDV  da  dn  ^ 
dx  dy  da  dt 


(3.9) 


Note  that  co  is  the  vertieal  veloeity  eomponent  in  the  sigma  eoordinate  system.  The 
horizontal  momentum  equations  are  given  by 
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(3.10) 
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dVD  dUVD  dV^D  dVa 
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(3.11) 


The  heat  and  salt  equations  are  represented  by 

dTD  dTUD  dTVD  dTco 
-  + - + - +  - 
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(3.12) 

(3.13) 


dt  dx  dy  da  da 
The  eddy  viseosity  Km  and  thermal  diffusivity  Kh  are  expressed  by 

where  /  is  the  mixing  length,  q  is  the  turbulent  speed,  and  Sm  and  Sh  are  stability 
funetions  whieh  are  analytieal  derived,  algebraie  relation  funetionally  dependent  upon 
dU I dz  ,  dV I dz  ,  dp! dz  ,q,  and  /.  The  level-2  turbulent  elosure  seheme  is  used; 
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W  +  F, 
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dD  dn 
+  a  —  +  — 
dt  dt 


(3.16) 


The  transformation  to  the  Cartesian  vertieal  veloeity  is  given  by 

jKi  rrf  ^ 

W  =  o)  +  U  a  —  -I - +V  a  —  +  - 

y  dx  dx)  [  dy  dy  J 

whieh  satisfies  the  no-flow  boundary  eonditions  at  the  surfaee  and  the  bottom. 

fn  (3.10)  and  (3.11),  should  be  subtraeted  from  p  to  form  p'  before  the 

integration,  is  generally  the  initial  density  field  whieh  is  area  averaged  on  z-levels 
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and  then  transferred  to  sigma  eoordinates  in  the  exaet  same  way  as  the  initial  density 
field.  This  proeedure  should  reduee  the  truneation  errors  assoeiated  with  the  ealeulation 
of  the  pressure  gradient  term  in  sigma  eoordinate  over  steep  topography  (see  Mellor  et 
ah,  1994  and  Mellor  et  al.  1998  for  evaluation  of  this  error  in  POM). 

POM  generally  uses  the  Smagorinsky  diffusivity  for  horizontal  diffusion  although 
a  eonstant  or  biharmonie  diffusion  ean  and  has  been  used  instead.  The  Smagorinsky 
formula  is 


A^=CAxAy^\vV  +  {VVf 


(3.17) 


where 


VV  +  (VVy  =  [dujdxf  +{^dvldx  +  duldy'f  +{^dvldy'f 


1/2 


Values  of  C  (the  HORCON  parameter)  in  the  range,  0.10  to  0.20  seem  to  work  well,  but, 
if  the  grid  spaeing  is  small  enough  (Oey  et  al,  1985a,  b),  C  ean  be  nil.  An  advantage  of 
the  Smagorinsky  relation  is  that  C  is  non-dimensional;  related  advantages  are  that 
deereases  as  resolution  improves  and  that  A^^  is  small  if  veloeity  gradients  are  small. 

The  equations,  governing  the  dynamies  of  eoastal  eireulation,  eontain  fast  moving 
external  gravity  waves  and  slow  moving  internal  gravity  waves.  It  is  desirable  in  terms  of 
eomputer  eeonomy  to  separate  the  vertieally  integrated  equations  (external  mode)  from 
the  vertieal  strueture  equations  (internal  mode).  This  teehnique,  known  as  mode  splitting 
(Simons,  1974;  Madala  and  Piaesek,  1977)  permits  the  ealeulation  of  the  free  surfaee 
elevation  with  little  saerifiee  in  eomputational  time  by  solving  the  veloeity  transport 
separately  from  the  three-dimensional  ealeulation  of  the  veloeity  and  the  thermodynamie 
properties. 

The  veloeity  external  mode  equations  are  obtained  by  integrating  the  internal 
mode  equations  over  the  depth,  thereby  eliminating  all  vertieal  strueture.  Thus,  by 
integrating  Equation  (3.9)  from  cr  =  -1  to  cr  =  0  and  using  the  boundary  eonditions,  an 
equation  for  the  surfaee  elevation  ean  be  written  as 
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(3.18) 


Sn  SUD  dVD  ^ 

—  + - + - =  0 

dt  Sx  Sy 


The  Courant-Friedrichs-Levy  (CFL)  computational  stability  condition  on  the 
vertically  integrated,  external  mode,  transport  equations  limits  the  time  step  according  to 


1 

1 

Sx^ 

-1/2 


(3.19) 


1  /9 

whereC,  =2(g//)  +U^^;  the  expected,  maximum  velocity.  There  are  other 

restrictions  but  in  practice  the  CFL  limit  is  the  most  stringent.  The  model  time  step  is 
usually  90%  of  this  limit.  The  internal  mode  has  a  much  less  stringent  time  step  since  the 
fast  moving  external  mode  effects  have  been  removed.  For  typical  coastal  ocean 
conditions  the  ratio  of  the  time  steps,  ,  is  often  a  factor  of  30-80  or  larger.  For 

more  information  on  the  sensitivity  of  POM  to  time  steps  (Ezer  et  ah,  2002).  In  SWAFS, 
the  horizontal  resolution  is  2  km  and  the  time  step  is  50  s. 

Together  with  the  turbulence  sub-model,  the  model  produces  realistic  bottom 
boundary  layers  which  are  important  in  coastal  waters  (Mellor,  1985)  and  in  tidally 
driven  estuaries  (Oey  et  ah,  1985a,  b)  which  the  model  can  simulate  since  it  does  have  a 
free  surface. 


SWAFS  is  run  using  a  24  hour  hindcast  and  48  hour  forecast  (Figure  11).  The 
hindcast  is  utilized  to  incorporate  ‘real-time’  data.  Unlike  atmospheric  models, 
observations  come  in  slower  and  are  not  as  spatially  diversified.  There  are  areas  within 
the  ocean  that  very  few  observations  of  the  subsurface  are  ever  recorded. 
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SWAPS  Assimilation  of  Obsei'v  ed  and  Synthetic 
T,S  Profiles  for  3  day  run 


Same  Obs  but  different  time  correlations 


Figure  11.  SWAPS  model  run  eyele.  (From;  Haeger,  S.D.,  2006) 

The  data  incorporated  into  the  hindcast  originates  from  a  variety  of  sources.  These 
sources  can  be  generalized  as  falling  into  one  of  the  following  categories:  (1)  sub-surface 
assimilation,  (2)  near-surface  assimilation,  (3)  lateral  boundary  conditions,  (4) 
meteorological  fluxes,  and  (5)  tides.  Each  of  these  categories  has  a  different  path  to  get 
data  into  SWAPS.  Data  assimilation  for  sub-surface  assimilation  starts  with  near  real¬ 
time  altimetry.  Sea  Surface  Temperature  (SST)  and  a  database  of  climatological  data 
being  ingested  into  the  Modular  Ocean  Data  Assimilation  System  (MODAS).  MODAS 
then  creates  synthetic  Conductivity/Temperature/Depth  (CTD)  profiles,  which  are  then 
added  to  the  spotty  real  CTD’s  and  Bathythermograph  (BT)  data.  These  profiles  are  then 
used  to  create  temperature  profiles  and,  indirectly,  salinity  profiles  through  optimum 
interpolation  (01).  This  data  is  then  spatially  ‘relaxed’,  or  decorrelated,  to  prevent 
‘bullseyes’  from  appearing  in  the  output.  Near-surface  data  assimilation  has  a  much 
shorter  path,  which  primarily  consists  of  other  multichannel  SSTs  (MCSSTs)  that  are 
incorporated  in  via  01.  This  data  is  decorrelated  in  much  the  same  way  as  discussed 
previously.  Tides  can  be  utilized  by  enabling  coastal  stations.  This  is  done  if  it  is 
determined  that  the  region  benefits  from  tidal  information.  This  is  not  the  case  for  all 

regions,  especially  the  open  ocean.  Meteorological  Fluxes  are  incorporated  by  taking  a 
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NOGAPS-COAMPS  blend  into  the  model.  As  SWAPS  goes  through  a  daily  run,  a  24 
hour  hindcast  is  run  with  the  previously  mentioned  data  getting  assimilated.  As  the 
forecast  portion  of  the  model  run  begins,  the  same  data  will  be  incorporated  in  again  with 
less  spatial  dependence.  As  seen  in  Figure  11,  this  is  done  again  at  the  24  hour  forecast 
mark.  The  reason  this  is  done  is  that  through  experimentation  by  the  modelers,  this 
method  has  been  determined  to  give  the  best  results.  Several  of  the  features  discussed 
previously  can  be  utilized  or  turned  off  depending  on  the  region  that  SWAPS  is  being 
run.  In  the  confidence  level  assessment  of  MOD  AS  (CLAM)  study,  models  are  tested  in  a 
given  region  to  determine  which  features  work  best  (Fox,  2000).  Once  the  best 
configuration  has  been  determined,  SWAPS  will  be  run  as  discussed  previously. 

C.  HPAC  AND  SECOND-ORDER  CLOSURE  INTEGRATED  PUFF 

(SCIPUFF) 

The  following  overview  is  taken  directly  from  Sykes  et  al  (2000).  It  is  used  to 
accurately  predict  the  effects  of  hazardous  material  releases  into  the  atmosphere  and  its 
impact  on  civilian  and  military  populations.  The  system  uses  integrated  source  terms, 
high-resolution  weather  forecasts  and  particulate  transport  analyses  to  model  hazard  areas 
produced  by  military  or  terrorist  incidents  and  industrial  accidents.  The  HPAC  system 
can  also  help  answer  the  question,  “How  good  is  the  prediction?”  by  providing 
probabilistic  calculations.  The  hazard  area  feature  estimates  the  weather  uncertainty  and 
turbulence  effects  on  possible  plume  trajectories  and  calculates  the  areas  of  hazard  impact 
and  the  degree  of  confidence  of  the  prediction. 

This  is  where  the  HPAC  model  is  unique  because  it  can  provide  a  hazard 
prediction  with  a  corresponding  estimate  of  confidence.  It  bases  this  confidence  on 
estimates  of  the  uncertainty  inherent  in  the  forecast  or  observed  weather  data.  These 
estimates  are  calculated  using  real  time  probabilistic  methods  or  via  empirical  models 
embedded  within  the  software.  The  SCIPUFF  model  is  the  atmospheric  transport 
empirical  model  used  in  HPAC. 

Transport  of  contaminants  in  the  turbulent  atmosphere  has  an  intrinsic  stochastic 
character.  Typical  deterministic  predictive  models  provide  an  ensemble  mean 
concentration  value,  but  are  not  able  to  predict  the  exact  path  and  diffusion  due  to 
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unpredictable  turbulent  motions  of  the  atmosphere.  A  probabilistic  transport  model 
predicts  both  ensemble  means  and  variances,  and  has  a  unique  role  in  that  it  can  provide 
rational  uncertainty  bounds  on  prediction  of  contaminant  concentrations  and  dosages. 

SCIPUFF  model  is  the  atmospheric  transport  engine  currently  used  in  HP  AC. 
SCIPUFF  is  a  Lagrangian  transport  and  diffusion  model  for  atmospheric  dispersion 
applications.  SCIPUFF  describes  two  basic  aspects  of  the  model.  First,  the  numerical 
technique  employed  to  solve  the  dispersion  model  equations  is  the  Gaussian  puff  method 
(Bass,  1980)  in  which  a  collection  of  three-dimensional  puffs  is  used  to  represent  an 
arbitrary  time-dependent  concentration  field.  Second,  the  turbulent  diffusion 
parameterization  used  in  SCIPUFF  is  based  on  the  second-order  turbulence  closure 
theories  of  Donaldson  (1973)  and  Lewellen  (1977),  providing  a  direct  connection 
between  measurable  velocity  statistics  and  the  predicted  dispersion  rates. 

The  Lagrangian  puff  methodology  affords  a  number  of  advantages  for 
atmospheric  dispersion  applications  from  localized  sources.  The  Lagrangian  scheme 
avoids  the  artificial  diffusion  problems  inherent  in  any  Eulerian  advection  scheme,  and 
allows  an  accurate  treatment  of  the  wide  range  of  length  scales  as  a  plume  or  cloud  grows 
from  a  small  source  size  and  spreads  onto  larger  atmospheric  scales.  This  range  may 
extend  from  a  few  meters  up  to  continental  or  global  scales  of  thousands  of  kilometers.  In 
addition,  the  puff  method  provides  a  very  robust  prediction  under  coarse  resolution 
conditions,  giving  a  flexible  model  for  rapid  assessment  when  detailed  results  are  not 
required.  The  model  is  highly  efficient  for  multiscale  dispersion  problems,  since  puffs 
can  be  merged  as  they  grow  and  resolution  is  therefore  adapted  to  each  stage  of  the 
diffusion  process. 

The  efficiency  of  SCIPUFF  has  been  improved  by  the  implementation  of  adaptive 
time  stepping  and  output  grids.  Each  puff  uses  a  time  step  appropriate  for  resolving  its 
local  evolution  rate,  so  that  the  multiscale  range  can  be  accurately  described  in  the  time 
domain  without  using  a  small  step  for  the  entire  calculation.  The  output  spatial  fields  are 
also  computed  on  an  adaptive  grid,  avoiding  the  need  for  the  user  to  specify  grid 
information  and  providing  a  complete  description  of  the  concentration  field  within  the 

computational  constraints  under  most  conditions. 
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The  generality  of  the  turbulenee  elosure  relations  provides  a  dispersion 
representation  for  arbitrary  eonditions.  Empirieal  models  based  on  speeifre  dispersion 
data  are  limited  in  their  range  of  applieation,  but  the  fundamental  relationship  between 
the  turbulent  diffusion  and  the  veloeity  fluetuation  statisties  is  applieable  for  a  mueh 
wider  range.  Understanding  of  the  daytime  planetary  boundary  layer  veloeity  fluetuations 
provides  reliable  input  for  the  second-order  closure  description  of  dispersion  for  these 
conditions.  The  closure  model  has  been  applied  on  local  scales  up  to  50  km  range  (Sykes 
et  al.,  1988)  and  also  on  continental  scales  up  to  3000  km  range  (Sykes  et  al.,  1993). 


The  second-order  closure  model  also  provides  the  probabilistic  feature  of 
SCIPUFF  through  the  prediction  of  the  concentration  fluctuation  variance.  In  addition  to 
giving  a  mean  value  for  the  concentration  field,  SCIPUFF  provides  a  quantitative  value 
for  the  random  variation  in  the  concentration  value  due  to  the  stochastic  nature  of  the 
turbulent  diffusion  process.  This  uncertainty  estimate  is  used  to  provide  a  probabilistic 
description  of  the  dispersion  result,  and  gives  a  quantitative  characterization  of  the 
reliability  of  the  prediction.  For  many  dispersion  calculations,  the  prediction  is  inherently 
uncertain  due  to  a  lack  of  detailed  knowledge  of  the  wind  field  and  a  probabilistic 
description  is  the  only  meaningful  approach. 

1.  Gaussian  Moment  Definition 

SCIPUFF  uses  a  Gaussian  puff  representation  for  the  concentration  field  of  a 
dispersing  contaminant.  A  three-dimensional  Gaussian  is  completely  described  by  its 
spatial  integral  moments  up  to  second-order,  and  can  be  written  in  the  form 


c(x)  = 


Q 


{27iy^{Det{a)y 


-17  exp 


(3.21) 


For  atmospheric  dispersion  problems,  one  must  also  consider  the  effects  of  the  ground 
surface  and  the  capping  inversion  at  the  top  of  the  planetary  boundary  layer,  which  are 
usually  represented  as  reflective  surfaces.  The  specific  Gaussian  variation  (3.21)  applies 
to  an  individual  puff,  but  in  general  the  local  concentration  field  will  be  composed  of  a 
sum  of  contributions  from  a  number  of  such  puffs. 


30 


2,  Moment  Transport  Equations 

The  advection-diffusion  equation  for  a  scalar  quantity  in  an  incompressible  flow 
field  can  be  written  as 

—  +—iuf)  =  kW^c  +  S  (3.22) 

dt  Sx. 

where  is  the  turbulent  velocity  field,  A:  is  the  molecular  diffusivity,  and  S 

represents  the  source  terms. 

The  atmospheric  velocity  field  is  generally  turbulent;  therefore  the  Reynolds 
averaging  technique  is  used  to  define  a  mean  and  a  turbulent  fluctuation  value.  Denoting 
the  mean  by  an  overbar  and  the  fluctuation  by  a  prime,  one  getsw  =u  +  u' ,  and  a  similar 
decomposition  can  be  applied  to  the  concentration,  c.  The  Reynolds  averaged 
conservation  equation  for  the  mean  scalar  concentration  is  thus 

—  +  —  (u.c)  =  -  —  (uic')  +  kV^c  +  S  (3.23) 

St  sx,^  >  5 xy  > 

where  -u'.c' is  the  turbulent  concentration  flux. 

3.  Turbulence  Closure  Diffusion  Model 

The  diffusion  model  in  SCIPUFF  is  based  on  second-order  turbulence  closure, 
which  provides  a  transport  equation  for  the  second  order  fluctuation  terms.  First-order 
closure  prescribes  the  turbulent  fluxes  in  terms  of  the  local  mean  gradients  using  an 
empirical  turbulent  diffusivity,  but  a  more  general  relation  can  be  obtained  from  a  higher- 
order  closure.  It  is  not  the  intention  to  review  turbulence  closure  theory  here,  only 
provide  the  basic  model  description  is  provided;  the  interested  reader  can  find  detailed 
discussions  in  the  literature,  e.g.,  Mellor  and  Herring  (1973),  Launder  et  al.  (1975), 
Lewellen  (1977). 
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The  equations  for  the  puff  moments  involve  the  turbulent  flux  of  eoneentration, 
and  a  rigorous  eonservation  equation  ean  be  derived  for  this  quantity  from  the  sealar  and 
momentum  equations.  Negleeting  the  moleeular  diffusion  terms,  the  flux  transport 
equation  ean  be  written  in  the  form 


s—  -  S 

- M,  C  +U. - 

St  ^  5x, 


Uf  C  =U-U, 
‘  ‘  J 


I  Sc  I  ,  Su^ 


Sx 


--U,  c 

Sx. 


—  u'u'c'-c'^  +  ^c'0'  (3.24) 
Sx,  ‘  ^  Sx,  n 


This  equation  involves  higher-order  terms,  sueh  as  the  triple  eorrelation  and  the 
pressure  eorrelation,  whieh  must  be  modeled  empirieally.  The  model  of  Lewellen  (1977) 
is  used  to  form  a  elosed  equation  for  the  turbulent  flux,  giving 
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(3.25) 


where  A  =  0.75  and  0.3  are  empirieal  model  eonstants.  The  turbulent  veloeity  seale, 

q,  is  defined  asq^  =u'u' ,  and  is  the  turbulent  length  seale;  these  quantities  will  be 
diseussed  more  fully  below.  The  buoyaney  term  uses  the  Boussinesq  approximation,  and 
g.  =(0,0,g)is  the  gravitational  aeeeleration,  T^is  the  referenee  temperature,  and  0'  is 

the  potential  temperature  fluetuation. 

The  general  equation  (3.25)  for  the  turbulent  fluxes  ean  be  integrated  spatially  to 
provide  transport  equations  for  the  flux  moments.  However,  the  multiple  tensor  indiees 
give  a  large  number  of  eorrelations  to  eonsider  in  the  general  ease.  Therefore,  some 
restrietions  are  introdueed  for  the  atmospherie  dispersion  oases  that  allow  this  study  to 
negleot  many  of  the  eorrelations.  First,  only  eonsider  the  vertioal  oomponent  of  the 

turbulent  drift,  since  advection  by  the  mean  wind  will  generally  dominate  in  the 

horizontal  directions.  Second,  the  only  off-diagonal  component  of  the  flux  moment 

tensor,  ,  to  be  considered  will  be  the  symmetric  horizontal  term; 
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X  U-,  c 


+ 


x,  u.  c 


(3.26) 


In  general,  the  off-diagonal  terms  represent  puff  distortions  due  to  veloeity  eovarianees. 

4.  Concentration  Fluctuation  Variance 

The  probabilistie  aspeet  of  the  SCIPUFF  dispersion  predietion  is  based  on  the 
seeond-order  elosure  model  for  the  eoneentration  fluetuation  varianee.  The  dispersion  of 
any  speeies  in  a  turbulent  veloeity  field  is  a  random  proeess  sinee  the  turbulent 
fiuetuations  are  effeetively  ehaotie  and  eannot  be  measured  or  predieted  in  detail.  The 
sealar  eoneentration  is  therefore  a  stoehastie  quantity,  with  a  probability  distribution  that 
depends  on  the  distribution  of  veloeity  fiuetuations.  Traditional  deterministie  estimates  of 
atmospherie  dispersion  only  provide  a  single  eoneentration  value  as  a  funetion  of  spaee 
and  time,  and  this  eorresponds  to  the  mean  value,  c,  for  some  definition  of  the  statistieal 
ensemble.  The  mean  value  is  the  first  moment  of  the  probability  distribution,  and 
eontains  no  information  about  the  statistieal  variability  in  the  predietion.  Higher  moments 
are  required  to  give  a  quantitative  deseription  of  the  variability.  The  probabilistie  aspeet 
of  the  SCIPUFF  dispersion  predietion  is  based  on  a  transport  equation  for  the  statistieal 
varianee  in  the  eoneentration  value  that  is  the  seeond  moment  of  the  probability 
distribution.  The  mean  and  the  varianee  are  then  used  to  provide  a  probabilistie  predietion 
using  a  parameterized  probability  density  funetion. 

The  key  aspeet  of  the  sealar  varianee  predietion  is  the  dissipation  timeseale. 
Proper  eharaeterization  of  the  dissipation  timeseale  allows  an  aeeurate  predietion  of  the 
eoneentration  fluetuation  varianee,  as  has  been  demonstrated  in  eomparison  with 
laboratory  data  (Sykes,  Lewellen,  and  Parker,  1986)  and  also  with  large-seale 
atmospherie  dispersion  observations  (Sykes  et  ah,  1993).  The  seeond-order  elosure  model 
for  the  dissipation  rate  was  originally  developed  by  Sykes  et  al.  (1984)  using  the 
laboratory  data  of  Faekrell  and  Robins  (1982).  A  fundamental  diseovery  in  that  study  was 
that  the  sealar  dissipation  seales  were  an  internal  property  of  the  sealar  field  itself,  rather 
than  being  determined  exelusively  by  the  veloeity  fiuetuations.  Essentially,  veloeity 
eddies  with  seales  larger  than  the  sealar  plume  or  eloud  will  meander  the  entire  sealar 
distribution  but  not  eause  any  nonlinear  easeade  of  sealar  fluetuation  varianee  onto 
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smaller  scales.  The  turbulent  cascade  process  is  only  driven  by  eddies  with  scales  similar 
to  the  instantaneous  plume  size.  This  is  not  the  case  for  persistent  mean  wind  shear, 
which  can  distort  a  cloud  by  means  of  the  continued  stretching. 

The  scalar  variance  dissipation  model  introduced  by  Sykes  et  al.  (1984,  1986) 
uses  a  single  estimate  of  the  velocity  and  length  scales  to  define  the  dissipation.  This  is 
appropriate  for  the  case  of  plume  dispersion  in  two  dimensions  with  near-isotropic 
turbulence.  Under  the  more  general  conditions  of  larger  scale  horizontal  dispersion  and 
arbitrary  sources,  characterization  of  the  different  dissipation  rates  associated  with  each 
direction  are  required.  The  vertical  direction  must  be  distinguished;  two  horizontal  scales 
are  also  needed.  The  horizontal  information  is  needed  to  distinguish  between  plume -type 
sources,  where  the  diffusion  occurs  in  two  spatial  dimensions,  and  puff-type  sources, 
which  diffuse  in  all  three  dimensions.  As  part  of  the  concentration  fluctuation  prediction, 
therefore  three  length  scales  in  addition  to  the  fluctuation  variance  are  required. 

5,  Numerical  Techniques 

The  puff  moment  evolution  equations  given  in  Section  C.l  generally  increase  the 
size  of  the  puff  through  turbulent  diffusion  and  elongation  along  the  direction  of  the  wind 
shear.  As  the  puff  grows,  the  local  representation  of  the  turbulence  and  velocity  fields 
using  the  puff  centroid  location  becomes  increasingly  inaccurate.  When  the 
meteorological  fields  are  inhomogeneous,  the  accuracy  of  the  calculation  can  only  be 
maintained  by  splitting  puffs  into  smaller  components  that  can  sample  the  variations  in 
the  meteorology  explicitly.  A  grid-based  method  for  splitting  puffs  was  presented  by 
Sykes  and  Henn  (1992),  where  the  moment  method  of  Egan  and  Mahoney  (1972)  was 
extended  to  include  some  shear  effects.  In  this  scheme,  grid  cells  can  be  thought  of  as 
each  containing  a  Lagrangian  puff.  As  the  puffs  move  and  spread  into  neighboring  grid 
cells  in  a  time  step,  the  masses  are  redistributed  so  as  to  maintain  a  single  puff  within 
each  cell.  To  avoid  any  numerical  grid  in  the  current  method  the  redistribution  cannot  be 
based  on  the  rectangular  grid  cells. 

The  three  diagonal  moments,  ,  give  the  effective  length  of  the  puff  along  each 
coordinate  direction,  and  then  the  puff  is  split  in  the  x-direction,  for  example,  when 
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CTjj  >  .  Here  Mi  represents  the  limit  for  the  horizontal  spread  and  should  be  chosen 

so  that  the  linearization  of  the  velocity  and  turbulence  fields  is  valid  for  a  Gaussian  with 
smaller  spread.  Similar  splitting  criteria  apply  for  the  other  two  coordinate  directions, 
using  MJ  in  the  y-direction  and  AV  in  the  z-direction. 

To  represent  the  original  Gaussian  puff  with  split  in  the  x-direction  as  several 
smaller,  overlapping  puffs  that  conserve  all  the  puff  moments  and  only  change  local 
concentration  values  by  a  small  amount,  the  original  puff  is  replaced  by  two  smaller  puffs 
in  the  following  manner.  The  new  centroid  locations  are  displaced  by  a  fraction,  r,  of  the 
puff  spread  in  the  x-direction,  and  by  a  distance  proportional  to  the  off-diagonal  moment 
in  the  other  two  coordinate  directions.  Thus 


=  xi  ±r 


(3.27) 


=  X2  ±r 


(3.28) 


-(«) 

X3 


=  X3  ±r 


(3.29) 


where  « =  {1,2}  corresponds  to  the  plus  and  minus  sign,  respectively.  The  diagonal 
moments  for  the  new  puffs  are  obtained  using  the  following  relations 
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and  the  new  off-diagonal  moments  are 

(3-33) 

^1?  =^13(1-''^)  (3-34) 
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(3.35) 


(a) 

^23 

Conservation  of  all  puff  moments  is  ensured  by  this  proeedure,  and  the  moments 
of  the  new  puffs  are  all  reduced  and  are  also  realizable.  The  realizability  constraints 
involve  the  Schwartz  inequality  between  the  diagonal  and  off-diagonal  moments.  These 
properties  can  be  verified  directly  from  the  above  relations. 

The  generalized  puff  description  contains  other  information  in  addition  to  the 
moments  discussed  above.  The  additional  variables  fall  into  two  categories,  either  a 
conserved  puff  integral  property  (similar  to  the  puff  mass)  or  a  puff  value  property  (such 
as  turbulence  length  scale).  Integral  properties  are  simply  divided  equally  between  the 
two  new  puffs,  and  the  value  properties  are  assigned  equally  to  both. 

6,  Meteorology  Specification 

To  make  accurate  calculations  and  predictions,  hazard  assessment  planning  and 
response  tools  must  have  access  to  three  general  types  of  weather  data:  (1)  historical 
weather  (climatology);  (2)  forecast  weather  (numerical  weather  predictions);  and  (3) 
current  weather  (observations).  The  type  of  weather  data  required  depends  upon  the 
needs  of  the  user.  To  plan  scenarios  that  are  either  notional  or  seven  days  previous  to  an 
examined  event,  historical  weather  or  climatology  data  is  most  suitable. 

To  plan  or  anticipate  events  within  a  seven-day  period,  forecast  weather  products 
are  the  data  of  choice.  For  a  post-event  or  real-time  response,  the  user  usually  requires 
access  to  current  observations.  Forecast  and  current  weather  data  is  available  to  the 

Hazard  Prediction  and  Assessment  Capability  program  via  the  Defense  Threat  Reduction 
Agency  (DTRA)  Meteorological  Data  Servers  (MDS),  while  historical  data  is  included  as 
part  of  the  tool  installation. 

The  meteorological  data  utilized  for  this  research  was  twofold.  For  the 
climatological  part  of  the  scenario,  the  climatology  inherent  in  HPAC  was  used.  For  the 
“real-time”  scenarios,  data  collected  from  COAMPS  forecasts  was  used. 
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D,  CHEMMAP^'^ 

The  following  overview  is  taken  direetly  from  MeCay  (2001).  Chemieal  spill 
modeling  provides  a  powerful  quantitative  tool  for  estimating  fate  and  potential  impaets 
of  ehemieal  releases.  Applied  Seience  Associates  (ASA)  has  been  developed  a  chemical 
spill  model,  CHEMMAP™,  to  predict  the  trajectory  and  fate  of  a  wide  variety  of 
chemical  products,  including  floating,  sinking,  soluble  and  insoluble  chemicals  and 
product  mixtures.  CHEMMAP™  simulated  a  number  of  processes  including:  (1)  slick 
spreading,  transport,  and  entrainment  of  floating  materials,  (2)  transport  of  dissolved  and 
particulate  materials  in  three  dimensions,  (3)  evaporation  and  volatilization,  (4) 
dissolution  and  adsorption,  (5)  sedimentation  and  resuspension,  (6)  and  degradation.  The 
model  uses  physical-chemical  properties  to  predict  the  fate  of  a  chemical  spill.  These 
include  density,  vapor  pressure,  water  solubility,  environmental  degradation  rates, 
adsorbed/dissolved  partitioning  coefficients  ),  viscosity,  and  surface  tension. 

1.  CHEMMAP™  Chemical  Database 

The  CHEMMAP^'^  chemical  database  contains  all  the  physical-chemical 
properties  required  by  the  model  to  simulate  the  transport  and  fate  of  the  spilled  material. 
The  database  includes  a  variety  of  text  and  numeric  descriptors  (names,  synonyms, 
registry  numbers,  etc.),  with  which  the  user  may  identify  and  characterize  the  spilled 
chemical  or  mixture.  The  model  is  capable  of  simulating  spills  of  pure  chemicals, 
chemicals  in  aqueous  or  hydrophobic  solutions,  or  chemicals  in  emulsions  (i.e.,  mixtures 
of  particulate  material  suspended  in  an  aqueous  base).  Thus,  the  database  also  includes 
characteristics  that  define  these  mixtures  and  solutions.  CHEMMAP^'^  uses  either  the 
Chemical  Abstract  System  (CAS)  registry  number  or  the  ETN  number  to  index  the 
chemical.  In  addition,  two  other  codes  uniquely  define  the  chemical  state  and  formulation 

to  be  simulated:  (1)  State  under  spill  conditions,  (2)  Eormulation  number  (to  allow 
multiple  database  entries  for  different  concentrations  of  the  same  chemical  and  state 
code). 

The  chemical  database  is  contained  in  a  Microsoft  Access  table.  The  data  may  be 
viewed  and  edited  from  the  user  interface.  The  initial  database  of  approximately  470 
chemicals  was  taken  from  Trench  et  al.  (1996).  These  data  will  continue  to  be  updated  as 
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new  information  becomes  available.  At  present,  the  database  has  been  increased  to 
approximately  900  chemical  substances,  including  those  commonly  transported  in  bulk 
and  most  likely  to  be  spilled. 

Several  properties  vary  with  temperature.  Thus,  the  database  values  are  for  a 
standardized  temperature  of  25“C.  The  model  corrects  these  parameters  to  the  ambient 
temperature  for  the  spill  incident.  The  algorithms  for  changing  viscosity  and  vapor 
pressure  to  ambient  temperature  are  taken  from  French  et  al.  (1996).  For  pure  chemical 
processes,  the  increase  per  10°C  is  assumed  2. 

2.  Chemical  Fates  Model 

The  chemical  fates  model  estimates  the  distribution  of  chemical  (as  mass  and 
concentrations)  on  the  water  surface,  on  shorelines,  in  the  water  column  and  in  the 
sediments.  The  model  is  three-dimensional,  separately  tracking  surface  slicks,  entrained 
droplets  or  particles  of  pure  chemical,  chemical  adsorbed  to  suspended  particulates,  and 
dissolved  chemical.  Processes  that  are  simulated  include  spreading,  transport,  dispersion, 
evaporation-volatilization,  entrainment,  dissolution,  partitioning,  sedimentation,  and 
degradation. 

The  model  initializes  the  spilled  mass  at  the  location  and  depth  of  the  release,  in  a 
state  dependant  upon  the  physical-chemical  properties  of  the  material.  The  state  code  and 
solubility  are  the  primary  determining  factors  for  the  initialization  algorithm.  If  the 
chemical  is  highly  soluble  in  water,  and  is  either  a  pure  chemical  or  dissolved  in  water 
(before  it  is  spilled),  the  chemical  mass  is  initialized  in  the  water  column  in  the  dissolved 
state. 

If  the  chemical  is  an  insoluble  liquid  or  a  gas,  the  following  initialization  is  used. 
If  the  density  is  less  than  or  equal  to  that  of  water  and  the  release  is  at  or  above  the  water 
surface,  the  model  initializes  the  material  in  surface  slicks.  Otherwise,  the  spilled  mass  is 
initialized  in  the  water  column  at  the  release  depth  in  a  plume  volume  defined  as  for  the 
soluble  chemicals.  Once  dissolution  has  occurred,  chemical  in  the  water  column  is  treated 
as  partially  adsorbed  to  suspended  sediment  particles  and  partially  dissolved.  Partitioning 
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between  these  states  is  assumed  to  be  in  eonstant  proportions  (based  on  linear  equilibrium 
theory,  using  from  the  ehemieal  database). 

Solids,  whieh  are  spilled  in  a  partieulate  form,  dissolve  and  then  are  partitioned  in 
the  water  eolumn  between  dissolved  and  partieulate  using  equilibrium  partitioning 
theory.  This  proeess  allows  for  a  stoehastie  approaeh  to  a  Mine  Drift  Seenario. 

Chemieal  mass  is  transported  in  three-dimensional  spaee  and  time,  by  surfaee 
wind  drift,  other  eurrents,  and  vertieal  movement  in  aeeordanee  with  buoyaney  and 
dispersion.  The  model  simulates  adsorption  onto  suspended  sediment,  resulting  in 
sedimentation  of  material.  Stokes  Law  is  used  to  eompute  the  vertieal  veloeity  of  pure 
ehemieal  partieles  or  suspended  sediment  with  adsorbed  ehemieal. 

Wind-driven  eurrent  (drift)  is  ealeulated  within  the  fates  model,  based  on  hourly 
wind  speed  and  direetion  data.  Surfaee  wind  drift  of  oil  has  been  observed  in  the  field  to 
be  1-6%  of  wind  speed  in  direetions  within  0-30  degrees  to  the  right  (in  the  northern 
hemisphere)  of  the  down-wind  direetion  (Youssef  and  Spaulding,  1993). 

The  horizontal  turbulent  diffusion  (randomized  mixing)  eoeffieient  normally 
ranges  from  1-10  m^/see  in  coastal  and  marine  waters.  The  vertical  turbulent  diffusion 
(randomized  mixing)  coefficient  is  typically  0.0001-0.001  m^/sec .  These  are  reasonable 
values  based  on  Okubo  (1971). 

For  surface  slicks,  the  model  estimates  surface  spreading,  slick  transport, 
entrainment  into  the  water  column,  and  evaporation,  to  determine  trajectory  and  fate  at 
the  surface.  Spreading  is  simulated  using  the  algorithm  of  Fay  (1971).  Entrainment  is 
modeled  as  for  oil,  using  data  in  Delvigne  and  Sweeney  (1988).  Surface  slicks  interact 
with  shorelines,  depositing  and  releasing  material  according  to  whether  the  material  is 
sticky  and  to  shoreline  type.  The  algorithms  used  are  those  developed  for  oil  spills,  as 
described  in  French  et  al.  (1999). 

Volatilization  from  the  water  column  is  calculated  from  the  chemical’s  vapor 
pressure  and  is  a  strong  function  of  temperature.  Degradation  is  estimated  assuming  a 
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constant  rate  of  "decay"  specific  to  the  environment  where  the  mass  exists  (i.e., 
atmosphere,  water  column  or  sediment). 

The  fates  model  computes,  in  space  and  time,  the  following:  (1)  area  covered  by 
surface  slicks,  (2)  radius  and  thickness  of  surface  slicks,  (3)  total  concentration  in  the 
water  column,  (4)  dissolved  concentration  in  the  water  column,  (5)  area  and  length  of 
shorelines  contaminated,  (6)  mass  per  unit  area  on  the  shorelines.  These  model  output 
data  are  mapped  by  the  user  interface  in  an  animated  display. 

3.  Environmental  Data 

A  geographical  database  supplies  data  for  water  depth,  shoreline  type,  sediment 
type,  habitat  type,  and  ice  cover  throughout  the  domain.  The  geographical  area  around  the 
spill  site  is  gridded  using  the  shoreline  location  to  define  land  and  water  cells. 

The  user  supplies  an  hourly  wind  time  series  specific  to  the  time  and  location  of 
the  spill.  For  this  research,  10m  (surface)  winds  from  COAMPS  were  utilized.  This  data 
is  important  for  insoluble  floating  chemicals,  but  has  little  influence  on  dissolved 
contaminants. 

The  chemical  transport  model  is  designed  so  that  it  may  accept  current 
(hydrodynamic)  data  from  a  variety  of  sources.  For  this  research,  surface  currents  from 
SWAFS  were  utilized. 

4,  Stochastic  Model 

The  probability  of  impact  from  a  chemical  discharge  is  quantified  using  ASA’s 
three-dimensional  stochastic  model  in  CHEMMAP^^.  In  the  stochastic  model,  many  runs 
are  made  for  each  release  scenario  and  location.  Using  a  historical  wind  and  current 
record  for  a  location,  the  model  randomizes  spill  date,  and  thus  the  wind  and  current 
conditions  selected  for  the  run.  In  this  way,  historically  observed  environmental 
conditions  are  sampled.  The  spill  volume  may  also  be  randomized.  An  effects  endpoint  is 
selected  based  on  toxicity  data,  which  provides  a  threshold  above  which  probabilities  and 
expected  concentrations  are  reported. 
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Model  output  includes  maps  of  chemical;  (1)  mass  or  volume  on  the  water 
surface,  (2)  mass  or  volume  on  shorelines,  and  (3)  particulate  concentration  in  the  water, 
(4)  dissolved  concentration  in  the  water. 

Five  statistics  are  produced  for  each  location  (cell)  in  the  model  grid.  These  are 
displayed  in  contour  maps,  which  may  be  overlaid  on  maps  of  receptors  of  concern. 
Statistics  on  all  model  runs  may  be  viewed  in  total  (as  probabilities  or  maximum 
exposures)  or  individually.  The  outputs  are;  (1)  probability  of  (any)  amount  exceeding 
the  threshold  passing  that  location,  (2)  time  when  amount  first  exceeds  the  threshold  at 
the  location,  (3)  mean  expected  maximum  mass  or  concentration  at  the  location  (i.e., 
peak  exposure  in  time  and  for  mean  environmental  conditions),  (4)  worst-case  (maximum 
possible)  amount  that  could  hit  a  location  (i.e.,  peak  exposure  in  time  and  maximum  of  all 
runs  under  all  environmental  conditions),  and  (5)  start  date  and  time  for  run  producing 
the  worst  case  at  the  location. 

An  assessment  of  potential  operational  risks  may  be  based  on  these  modeling 
results.  The  mean  expected  and  worst-case  exposure  concentrations  can  be  compared  to 
effects  endpoints  of  interest.  The  worst-case  spill  for  any  given  location  may  be  evaluated 
in  more  detail.  The  stochastic  model  identifies  the  conditions  that  constitute  a  worst  case. 
That  scenario  may  then  be  run  with  the  (single  scenario)  chemical  fates  model  to  obtain 
animated  displays  of  concentrations  over  time.  These  fates  results  may  be  input  to  the 
biological  effects  model  to  quantitatively  estimate  impacts. 
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IV.  MATHEMATICAL  AND  STATISTICAL  ANALYSIS 


Before  one  ean  talk  about  analyzing  the  data,  the  methods  involved  to  analyze  the 
data  must  be  deseribed.  Both  outputs  of  SWAPS  and  COAMPS  had  to  be  evaluated 
separately.  Simple  statistieal  methods,  sueh  as  mean  and  varianee,  are  utilized  to 
eharaeterize  overall  properties  of  the  Northern  Persian  Gulf;  however,  to  get  an  idea  of 
the  general  flow  and  predominant  features  a  more  eomplex  method  had  to  be  used.  An 
Empirieal  Orthogonal  Funetion  (EOF)  method  is  used  to  find  temporal-spatial  variability. 
This  method  is  similar  to  Prineipal  Component  Analysis  (PC A)  in  statisties.  The 
following  seetions  will  deseribe  these  methods  and  will  break  down  how  these  methods 
are  ealeulated. 

A,  BASIC  MATHEMATICS  AND  STATISTICS 

One  of  the  most  basie  statistieal  methods  is  to  ealeulate  the  arithmetie  mean, 
whieh  is  often  simply  ealled  the  mean  or  average.  The  mean  is  the  sum  of  all  the  values 
divided  by  the  total  number  of  values,  as  in  the  following  equation: 

-  ^  (.^1  +X2  +X3  +  ...  +  x„  ^ 

n  n 


There  two  main  methods  of  using  the  mean  that  were  used  in  this  researeh.  One 
was  to  use  the  mean  to  ealeulate  the  perturbation,  or  ‘prime’  values.  This  is  done  by 
applying^  =  M -i-m' .  This  method  will  be  used  when  ealeulating  the  EOFs.  The  other 
method  was  used  to  get  a  general  idea  of  flows  in  the  Northern  Persian  Gulf  The  seale  of 
the  flow  and  general  eharaeter  of  the  flow  was  examined  using  this  method. 

In  addition  to  the  mean  and  average,  diseussed  early,  two  statistieal  methods  are 
utilized.  These  methods  inelude  standard  deviation  and  eovarianee.  Standard  deviation  of 
a  probability  distribution  is  defined  as  a  measure  of  the  spread  of  those  values.  It  is  also 
the  square  root  of  varianee.  Mathematieally,  standard  deviation  is  defined  as 


s  = 


(n-1) 


(4.2) 
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Standard  deviation  is  the  most  eommon  measure  of  statistical  dispersion, 
measuring  how  widely  spread  the  values  in  a  data  set.  If  the  data  points  are  close  to  the 
mean,  then  the  standard  deviation  is  close  to  zero.  As  the  data  points  move  farther  away 
from  the  mean,  the  standard  deviation  moves  farther  from  zero.  The  other  statistical 
method,  covariance,  is  the  measure  of  how  much  two  variables  vary  together.  This  is 
different  from  variance,  mentioned  early,  which  measures  how  much  a  single  variable 
varies. 


Covariance  is  computed  by 


cov  = 


i(x,-X)(Y,-Y) 


(n-1) 


(4.3) 


for  two  variables  X  and  Y.  If  two  variables  tend  to  vary  together  then  the  covariance 
between  the  two  variables  will  be  positive.  On  the  other  hand,  if  when  one  of  them  is 
above  its  expected  value,  the  other  variable  tends  to  be  below  its  expected  value,  then  the 
covariance  between  the  two  variables  will  be  negative.  If  the  variables  vary  independent 
from  each  other,  the  covariance  will  be  zero.  Conversely,  if  the  covariance  is  zero,  the 
variables  need  not  be  independent. 

Another  important  mathematical  method  involves  matrix  algebra.  Specifically  the 
covariance  matrix  and  eigenvectors/eigenvalues  of  a  given  matrix  will  be  discussed.  If 
the  data  set  has  more  than  2  dimensions,  there  is  more  than  one  covariance  measurement 
that  can  be  calculated.  A  useful  way  to  account  for  the  additional  variances  associated 
with  each  dimension  is  to  create  a  covariance  matrix.  This  matrix  is  composed  of  the 
covariances  between  elements  of  a  vector.  A  covariance  matrix  of  three  dimensions 
would  look  like 


^  cov(x,x) 

cov(x,y) 

cov(x,z)^ 

c  = 

cov(y,x) 

cov(y,y) 

cov(y,z) 

(4.4) 

l^cov(z,x) 

cov(z,y) 

cov(z,z) J 
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Some  points  to  note:  Down  the  main  diagonal,  the  eovarianee  value  is  between 
the  same  dimensions.  These  values  are  the  variances  for  the  given  dimension.  The  other 
point  is  that  sincecov(a,h)  =  cov(h,a) ,  the  matrix  is  symmetrical  about  the  main 
diagonal. 

Eigenvectors  are  a  special  case  of  multiplying  two  matrices  together.  There  are 
two  key  properties  of  all  eigenvectors.  First,  all  eigenvectors  multiplied  by  a  scalar  are 
equal,  since  all  the  scalar  does  is  make  the  vector  longer  but  not  change  its  direction. 
Secondly,  all  eigenvectors  of  a  matrix  are  perpendicular,  i.e.  at  right  angles  to  each  other, 
no  matter  how  many  dimensions  one  has.  Mathematically  another  term  for  perpendicular 
is  orthogonal.  This  is  important  because  it  means  the  data  can  be  expressed  in  terms  of 
these  perpendicular  eigenvectors,  instead  of  expressing  them  in  terms  of  the  x  and  y  axes. 
Another  important  thing  to  know  is  that  when  eigenvectors  are  found  it  is  common  to 
find  those  whose  length  is  exactly  one.  This  is  because  the  length  of  a  vector  does  not 
affect  whether  it  is  an  eigenvector  or  not,  whereas  direction  does.  So,  in  order  to  keep 
eigenvectors  standard,  whenever  an  eigenvector  is  found  it  is  scaled  to  have  a  length  of  1 , 
so  that  all  eigenvectors  have  the  same  length.  Eigenvalues  are  closely  related  to 
eigenvectors.  Eigenvectors  and  eigenvalues  always  come  in  pairs.  Eigenvalues  are  the 
multiple  used  to  get  the  eigenvector  back  to  the  original  length  it  was  before  being 
converted  to  a  length  of  one. 

B,  CONVENTIONAL  EOF  ANALYSIS 

Let  a  temporally  and  spatially  varying  scalar  variable  y/  be  represented  by 

with  {xi,  yj)  the  horizontal  grids,  Xk  the  time  sequence  in  months,  t,  = 
1,  2,  ...  ,  the  time  sequence  in  a  month.  Let  wi^i^yj^Zk)  t)e  the  monthly  mean  of  the 
scalar.  The  synoptic  anomalies 

¥  '(a  ^ypXk^^i)  =  wiZi  ,yj,Xk^'^i)-  wiXi  ,ypXk)^  (4-5) 

are  re-arranged  into  a  AxP  matrix,  ^  2,  ...,  A;  and  p  =  1,  2,  ...,  P.  Here  N 

is  the  total  number  of  the  horizontal  grid  points  and  P  is  the  total  number  of  time  points 
used  for  computing  the  covariance  matrix.  The  EOF  analysis  widely  used  in 
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oceanographic  and  meteorological  research  (e.g.,  Weare  et  ah,  1976;  Richman,  1986; 
Chu  et  ah,  1997  a,  b)  is  the  same  as  the  prineipal  eomponent  (PC)  analysis  (Hotelling, 
1933)  in  the  statistics  community.  PCs  are  the  amplitudes,  whieh  are  functions  of  time,  of 
their  eorresponding  EOFs.  These  EOFs  ean  be  found  by  calculating  the  unitary 
eigenvectors  of  the  covariance  matrix  associated  with  the  sample  data  field.  EOF  analysis 
separates  the  data  sets  into  eigenmodes.  Generally  speaking,  eaeh  mode  has  an  assoeiated 
variance,  dimensional  spatial  pattern,  and  non-dimensional  time  series.  The  NxN  spatial 
covariance  matrix  is  ealeulated  by 

7?ji  7?j2  ... 

_  ^21  ^22  •••  ^2N 

RnI  R-NI  •••  RnN 

where  n  and  m  (1,  2,  ...,  N)  denote  the  grid  locations.  The  diagonal  elements  of  the 
covariance  matrix  {Rnn}  are  the  varianee  at  loeation  r„.  The  off-diagonal  elements  are  the 

covariance  with  spatial  lag  equal  to  the  difference  between  the  row  and  column  indices. 
This  symmetric  matrix  has  N  real  eigenvalues /l^  ,  and  eigenveetors^zJ^(r^)  ,  such  that 

|;V<.(r)  =  4A(r,).  =  (4.7) 

7=1 

The  eigenvectors  ^zJj(r.),  •••’  are  called  EOFs.  Each  of  them  is  an  Ap¬ 

point  field  showing  anomaly  pattern.  The  eigenvalues,  (a  =1,2,...,  A^),  are  all  positive 
and  the  summation  of  them,  ,  equals  the  total  variance.  Therefore,  is  considered 

a 

as  the  portion  of  total  variance  ‘explained’  by  the  EOF  ^zJ„(r^)  .  It  is  convenient  to  label 
the  eigenfunctions  <I)„{y^)  so  that  the  eigenvalues  are  in  descending 

X^  >  X2  >  ...>  Xjy  (4.8) 

The  data  matrix,  ^'(r^,t^) ,  is  thus  approximately  written  by 


R.,n  =  X ^ ’  K >  )  ’  (4.6) 

^  ^  p=\ 
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(4.9) 


W  '(r„  ,tp)  =  Y^PC^  {tp  )(j)^  (r„ ) , 

a 

where  PC^it^)  is  the  principal  component  with  the  same  unit  as  the  scalar  y/  and  a  size 

of  P,  representing  the  temporal  variation  of  the  associated  spatial  pattern  described  by 
EOF^/i„(r„). 

C.  COMPLEX  EOF  ANALYSIS 

Consider  a  two-component  variable  such  as  the  horizontal  velocity  vector  (u,  v). 
A  complex  field  can  be  constructed  by 

w  =  u+iv,  /  =  (4.10) 

Similar  to  the  previous  section,  let  a  temporally  and  spatially  varying  complex 
field  w  be  represented  by  w{x.,yj,Zk,T,)  with  (x,,  yj)  the  horizontal  grids,  Zk  the  time 

sequence  in  months,  r,  =  1,  2,  ...,  n,  the  time  sequence  in  a  month.  Let  if(x.,y^.,;t'^)  be 
the  monthly  mean  of  the  scalar.  The  synoptic  anomalies 

w  '(a  ,yj,Z/c’h)  =  Mx,  ^yj^Xk)^  (4  •  1 1 ) 

are  re-arranged  into  a  NxP  matrix,  w'(r„,t^)  ,  n  =  1,  2,  ...,  A;  and  p  =  1,2,  ...,  P.  Here  N 

is  the  total  number  of  the  horizontal  grid  points  and  P  is  the  total  number  of  time  points 
used  for  computing  the  covariance  matrix.  Ax  A  spatial  covariance  matrix  is  calculated 
by 

i?ll  i?j2  ... 

_  ^21  ^22  •••  ^2N 

D  D  D 

_N\  ^N2  ^NN 

where  the  superscript  denotes  the  transpose.  The  diagonal  elements  of  the  covariance 
matrix  {Rnn}  are  the  variance  at  location  r„.  The  covariance  matrix  R  is  a  Hermitian,  so 
the  eigenvalues  { /l„  }  are  positive  real  numbers  whereas  the  EOFs  { }  are  complex, 

=  /  =  (4.13) 

7=1 


.  (4.12) 

^  ^  P=\ 
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Each  of  the  complex  EOFs^zJj(r.),  •••’  them  is  an  At-point  field 

showing  anomaly  velocity  vector  pattern.  The  eigenvalues,  A,^  {a  =\,2,...,N),  are  all 
positive  and  the  summation  of  them,  ,  equals  the  total  variance.  Therefore,  is 

a 

considered  as  the  portion  of  total  variance  ‘explained’  by  the  EOF  .  It  is  convenient 
to  label  the  eigenfunctions  so  that  the  eigenvalues  are  in  descending 

X,  >  X2  >  ...>  Xj^.  (4.14) 

The  data  matrix,  w'(r^,t^)  ,  is  thus  approximately  written  by 

^ '(r„ ,tp)  =  Y,raitp) exp[/6»„ (r„ )  ,  (4.1 5) 

a 

where  ]  are  the  temporally  varying  amplitude  (positive)  and  phase 

(between  -180°  and  180°)  associated  spatial  pattern  described  by  EOF  <I)^{y^)  . 
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V.  MINE  WARFARE  (MIW) 


Sea  mines  have  been  important  in  naval  warfare  throughout  history  and  continue 
to  be  so  today.  Originally  called  a  torpedo,  mines  have  caused  major  damage  to  naval 
forces,  slowed  or  stopped  naval  actions  and  commercial  shipping,  and  forced  the 
alteration  of  strategic  and  tactical  plans.  The  threat  posed  by  sea  mines  continues,  and  is 
increasing,  in  today’s  world  of  inexpensive  advanced  electronics,  nanotechnology,  and 
multiple  potential  enemies,  some  of  which  are  difficult  to  identify.  The  largely 
unregulated  sale  of  sea  mines  by  allies  and  third  parties  is  contributing  directly  to  this 
growing  threat. 

Because  of  the  low  cost  and  wide  availability  of  modern  sea  mines,  the 
importance  as  a  threat  to  shipping  and  naval  force  operations  is  growing  rapidly.  The 
threat  continues  to  increase  as  the  U.S.  naval  forces  operate  increasingly  in  the  world’s 
littorals.  The  Navy  is  responsible  for  protecting  all  maritime  forces,  including  logistics 
transport  and  Marine  Corps  units,  against  the  mine  threat  wherever  it  may  be 
encountered,  from  the  sea  lanes,  to  logistics  unloading  areas,  to  the  high-water  mark  on 
the  landing  beaches. 

Naval  mines  can  be  used  strategically,  channeling  or  denying  passage  through 
restricted  waters  and  in  and  out  of  ports  needed  for  sustenance  by  littoral  nations.  Mines 
can  shape  the  naval  battlespace,  the  approaches  to  it,  and  routes  of  commerce,  thus 
setting  the  conditions  of  a  campaign.  Used  tactically,  they  can  slow  or  stop  movement  to 
and  through  narrow  straits  and  to  landing  zones  on  beaches,  and  in  so  doing  can  also 
make  a  slowed  or  stopped  force  more  vulnerable. 

The  breakup  of  the  Soviet  Union,  with  its  enormous  stockpile  estimated  at  nearly 
half  a  million  sea  mines,  has  resulted  in  the  potential  widespread  availability  of  mines.  A 
vintage  World  War  I  mine  was  able  to  inflict  $96  million  worth  of  damage  to  the  USS 
Samuel  B.  Roberts  (FFG  58)  in  the  Persian  Gulf  (Zwolski,  1998).  There  were  several 
other  instances  in  the  late  1980’s  of  drifting  mines  being  deployed  around  shipping  lanes 
in  the  Persian  Gulf  region.  In  fact  75  percent  of  damage  to  U.S.  Navy  capitol  ships  from 
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1988  to  1998  came  from  mines,  two  of  whieh  were  World  War  I  teehnology  (Zwolski, 
1998).  Subsequently,  mines  have  gained  the  reputation  as  the  least  expensive,  yet  most 
effective,  offensive  and  defensive  weapon  of  war.  The  threat  is  still  there  today. 

In  the  eurrent  threat  environment,  a  mine  attaek  on  the  eeonomie  lifeline  of  Iraq  is 
anything  but  impossible.  Sinee  the  thesis  is  unclassified,  it  is  not  permitted  to  use 
elements,  which  would  be  used  in  the  case  of  an  aetual  attack;  however  using  a  simple 
drifter  can  give  a  good  impression  of  the  vulnerabilities  to  a  drifting  mine  attaek. 
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VI.  CHEMICAL  POLLUTANTS 


In  the  current  threat  environment,  a  chemical  attack  on  the  economic  lifeline  of 
Iraq  is  anything  but  impossible.  Since  the  thesis  is  unclassified,  the  chemicals  were 
chosen  due  to  abundance  in  the  region  and  do  not  necessarily  have  any  elements  which 
would  be  used  in  the  case  of  an  actual  weapon  with  mass  destruction  (WMD)  attack. 
These  chemicals  were  also  chosen  because  of  there  high  toxic  elements.  For  this  reason, 
the  choices  are  Ammonia  and  sulfur  mustard. 

A,  AMMONIA 

Ammonia,  NH3  or  NH4+,  is  prepared  commercially  in  vast  quantities,  with  major 
uses  of  Ammonia  and  its  compounds  being  fertilizer  and  refrigerant  and  in  the  production 
of  other  chemicals.  With  the  rise  in  natural  gas  prices,  fertilizer  production  has  shifted  to 
areas  with  cheaper  natural  gas  prices.  One  of  these  areas  is  the  Persian  Gulf  region.  As 
fertilizer  production  increase  in  the  Persian  Gulf  so  to  must  the  shipment  of  ammonia. 
Even  in  the  best  case  scenario,  increase  the  amount  that  is  shipped  increases  the  chances 
of  an  Ammonia  spill  within  the  Persian  Gulf  region. 

Ammonia  is  a  clear  colorless  gas  with  a  strong  odor.  Although  Ammonia  gas  is 
lighter  than  air,  vapors  from  a  leak  initially  hug  the  ground  (ASA’s  Database).  It  is 
shipped  as  a  liquid  under  its  own  vapor  pressure.  Its  density  in  the  liquid  form  is  12.8825 
kg/m  .  Contact  with  the  unconfined  liquid  can  cause  frostbite.  Gas  generally  regarded  as 
nonflammable  but  does  bum  within  certain  vapor  concentration  limits  and  with  strong 
ignition.  Fire  hazard  increases  in  the  presence  of  oil  or  other  combustible  materials. 

It  is  extremely  soluble  in  water;  one  volume  of  water  dissolves  about  1,200 
volumes  of  the  gas  at  0“C  (90  grams  of  Ammonia  in  100  cc  of  water),  but  only  about  700 
volumes  at  room  temperature  and  still  less  at  higher  temperatures  (Ammonia,  2007). 
Ammonia  can  be  tasted  in  water  at  levels  above  about  35  mg/liter.  Odor  thresholds  have 
been  variously  reported  according  to  the  definition  used  and  technique  of  measurement. 
Most  people  can  identify  Ammonia  in  air  at  about  35  mg/m  and  can  detect  it  at  about 
one-tenth  of  this  level. 
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Ammonia  is  a  floater,  highly  volatile,  highly  soluble  and  slightly  absorbable  to 
particles,  and  reacts  exothermically  with  all  acids.  Violent  reactions  are  possible.  It  also 
readily  combines  with  silver  oxide  or  mercury  to  form  compounds  that  explode  on 
contact  with  halogens.  When  in  contact  with  chlorates,  it  forms  explosive  ammonium. 

As  for  toxicity.  Ammonia  is  immediately  dangerous  to  life  or  health  indicator 
(IDLH)  is  relatively  small  (300  ppm),  its  short  term  exposure  limit  (STEL)  is  35  ppm  and 
its  odor  threshold  is  0.019  ppm.  Its  degradation  rate  is  0.1586  in  both  air  and  water. 
Contact  with  Ammonia  could  cause  skin  and  eye  burns  and  inhalation  some  burning 
sensation,  cough,  shortness  of  breath  and  sore  throat.  For  the  ecosystem,  its  slight  toxicity 
can  be  lethal  to  shrimp,  prawns,  salmon,  trout  and  catfish. 

B,  SULFUR  MUSTARD  /  MUSTARD  GAS 

Sulphur  mustards,  of  which  Mustard  Gas  is  a  member,  were  first  used  as  a 
chemical  warfare  agent  effectively  World  War  I  by  the  German  army  against  Canadian 
soldiers  in  1917  (Sidell,  et  al,  1998).  Mustard  Gas  is  chemically  called  beta-chloroethyl 
sulphide  (C4H8CI2S)  and  has  the  ability  to  form  large  blisters  on  exposed  skin.  The  affect 
of  blistering  puts  Mustard  Gas  in  a  subcategory  of  chemical  agents  called  blister  agents. 
This  blister  agent  is  very  popular  and  has  since  been  used  repeatedly  over  the  years.  This 
includes:  Iraq  against  Iran  in  1981  and  1983-1988,  Iraq  against  Iraq  in  1987-1988  and 
Iraq  against  Kurd  in  1988  (CBWInfo.com  Factsheets,  2007).  The  use  of  poison  gas, 
including  Mustard  Gas,  during  warfare  was  prohibited  by  the  Geneva  Protocol  of  1925 
and  the  subsequent  Chemical  Weapons  Convention  of  1993,  which  also  prohibited  the 
development,  production  and  stockpiling  of  such  weapons  (CBWInfo.com  Factsheets, 
2007). 

Despite  the  existence  of  more  highly  effective  chemical  warfare  agents  that  are 
more  toxic.  Mustard  Gas  has  not  lost  its  usefulness  because  of  its  special  characteristics. 
It  is  very  toxic  and  difficult  to  treat,  versatile,  persistent,  cheap,  easy  to  produce,  and 
difficult  to  protect  against.  Furthermore,  sulfur  mustard  is  toxic  as  droplets,  liquid  and 
vapor.  Due  to  the  simple  and  inexpensive  synthesis,  as  well  as  physical  properties,  it  has 
the  potential  of  being  one  of  the  leading  threats  from  both  terrorist  and  military  sources 
(Vidan,  2002). 
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Pure  sulfur  mustard  is  a  colorless,  viscous  liquid  with  a  slight  odor  of  castrol  oil  at 
room  temperature,  whereas  when  used  as  a  chemical  agent  it  is  typically  a  dark  yellow- 
brown  liquid  with  an  odor  typically  described  as  mustard,  horseradish  or  garlic  (Vidan, 
2002).  Sulfur  mustard  is  barely  soluble  in  water,  0.07%  at  10°C,  and  very  soluble  in 
organic  solvents,  fuels  and  lubricants  (Vidan,  2002).  In  aqueous  solutions,  mustard  agent 
decomposes  into  non-poisonous  products  by  means  of  hydrolysis.  This  reaction  is 
catalyzed  by  alkali.  However,  only  dissolved  mustard  agent  reacts,  which  means  that  the 
decomposition  proceeds  very  slowly. 

Sulfur  mustard  is  a  potent  alkylating  agent,  readily  penetrating  the  skin,  mucus 
membranes  and  the  eyes,  which  are  the  most  sensitive  (Vidan,  2002).  Most  victims 
sustain  non-fatal  but  disabling  injuries,  mostly  to  the  eyes,  respiratory  tract  and  skin. 
Ocular  injuries  appear  in  75±90%  of  all  Mustard  Gas  casualties,  with  reports  of  delayed 
ocular  morbidity  appearing  years  later  (Vidan,  2002).  Since  Mustard  Gas  primarily 
disables,  vice  kills,  for  every  individual  affected,  on  average,  another  must  be  utilized  to 
tend  to  the  affected  individual.  This  will  further  reduce  the  number  of  personnel  available 
to  complete  a  given  mission,  making  it  a  very  effective  chemical  agent. 

Since  density  of  sulfur  mustard  is  greater  than  water.  Therefore,  in  an  aquatic 
environment  the  primary  delivery  of  pure  sulfur  mustard  will  be  from  the  air.  To  use  it  as 
a  chemical  agent,  sulfur  mustard  will  be  disseminated  using  aerial  sprays  or  munitions 
and  will  be  initially  in  the  form  of  droplets.  These  droplets  would  be  subject  to 
gravitational  settling.  Sulfur  mustard  vapor  from  resulting  deposits  may  reenter  the 
atmosphere  by  evaporation.  In  the  atmosphere,  sulfur  mustard  vapor  will  degrade 
primarily  by  reaction  with  photochemically  produced  hydroxyl  radicals;  the  resulting  half 
life  being  estimated  as  1.4  days  (Spectrum  Fact  Sheet,  2007). 

When  sulfur  mustard  is  introduced  into  water  at  very  low  concentrations  so  that  it 
is  dissolved,  it  will  rapidly  hydrolyze.  At  25‘’C,  the  estimated  half-live  of  sulfur  mustard 
when  dissolved  in  large  amounts  of  water  are  4  minutes  (Spectrum  Fact  Sheet,  2007). 
Volatilization  will  be  slow  by  comparison.  Hydrolysis  in  seawater  will  be  a  factor  of  2.5 
slower  because  of  the  common-ion  effect  exerted  by  the  chloride  ion  (Spectrum  Fact 
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Sheet,  2007).  Despite  its  high  rate  of  hydrolysis,  sulfur  mustard  that  has  not  dissolved 
may  persist  for  longer  periods  of  time  beeause  its  rate  of  solution  is  slow. 

When  sulfur  mustard  is  introdueed  into  seawater  in  higher  eoneentrations,  three 
zones  with  different  behavior  are  observed  (Speetrum  Faet  Sheet,  2007).  A  small  amount 
of  sulfur  mustard  will  remain  on  the  surfaee  in  a  surfaee  film  where  it  will  be  destroyed 
by  hydrolysis  and  to  some  extent  volatilization;  this  surfaee  film  may  last  from  one  to 
several  days.  Bulk  sulfur  mustard,  or  droplets  formed  by  wind  aetion  on  the  surface  fdm, 
will  sink  to  the  bottom  of  the  water  because  sulfur  mustard  is  heavier  than  water.  This 
sulfur  mustard  will  be  degraded  by  hydrolysis;  however,  the  hydrolysis  rate  will 
essentially  be  the  rate  of  dissolution.  The  rate  of  dissolution  will  depend  on  the  exposed 
surface,  the  amount  of  agitation,  and  the  temperature.  In  the  water  column  of  a  body  of 
seawater,  dissolved  sulfur  mustard  will  be  lost  by  hydrolysis  with  its  half-life  being 
around  15  min  at  25°  C  (Spectrum  Fact  Sheet,  2007). 

Disposal  by  sea  burial  might  result  in  releases  of  large  quantities  of  sulfur 
mustard.  When  bulk  quantities  of  sulfur  mustard  are  disposed  of  at  sea,  some  of  it  will 
form  a  surface  film  that  will  disappear  as  a  result  of  hydrolysis  and  volatilization  within 
several  days  (Spectrum  Fact  Sheet,  2007).  Most  of  the  sulfur  mustard  will  sink  to  the 
bottom  of  the  water  because  it  is  heavier  than  water.  This  sulfur  mustard  will  be  degraded 
by  hydrolysis;  however,  the  hydrolysis  will  be  limited  by  the  rate  of  dissolution.  The  rate 
of  dissolution  will  depend  on  the  exposed  surface,  the  amount  of  agitation,  and  the 
temperature. 
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VII.  ANALYSIS  OF  PHYSICAL  MODEL  OUTPUT 


Datasets  were  eolleeted  via  two  different  operational  models,  COAMPS  and 
SWAPS.  The  following  discussion  will  examine:  (1)  where  the  data  came  from,  (2)  how 
it  was  collected,  (3)  what  methods  were  used  to  process  it,  both  spatially  and  temporally, 
and  (4)  how  the  two  datasets  compared  to  each  other  at  a  specific  geographical  reference. 

A,  PHYSICAL  MODEL  OUTPUT 

1.  Currents 

The  currents  dataset  was  retrieved  from  the  Operational  run  of  SWAPS.  This 
model  is  run  once  daily  at  Navy  Oceanographic  Office  (NAVO)  located  at  Stennis  Space 
Center  near  Bay  St.  Louis  Mississippi.  The  model  is  run  with  a  24  hour  hindcast  and  a  48 
hour  forecast.  Personnel  at  the  NAVO  removed  data  every  hour  from  the  first  day  of  the 
forecast  portion  of  the  model  and  stitched  these  Tau’s  together  into  month  long  datasets. 
The  data  for  currents  was  split  into  monthly  intervals  due  to  the  size  of  the  data  involved. 
On  average,  one  month  worth  of  currents  equaled  450  mb  of  data.  If  the  data  had  been 
collected  in  one  long  dataset  it  would  have  been  several  gigabytes  in  size.  This  in  turn 
would  have  created  issues  with  transferring  it  to  the  Naval  Postgraduate  School  and  in 
analyzing  the  data.  The  area  is  from  23.5°N  to  30.5°  N  and  47.5°E  to  57.5°E.  The 
horizontal  resolution  of  the  operational  model  was  2  km.  The  data  was  retrieved  in  hourly 
intervals  from  Eebruary  through  July  2006.  Each  hour  field  contained  both  U  and  V 
components  of  0.5  m  deep  currents,  referred  to  hereafter  as  currents.  There  were  a  few 
days  that  the  data  was  either  not  collected  or  was  erroneous.  This  data  was  removed 
during  the  data  analysis  process  and  replaced  with  the  mean  calculated  from  the  rest  of 
the  dataset. 

2,  Winds 

The  wind  dataset  was  retrieved  from  the  Operational  run  of  COAMPS.  This 
model  is  run  operationally  with  incremental  updates  throughout  each  day  from  Elect 
Numerical  Meteorology  and  Oceanography  Center  (ENMOC)  in  Monterey  California. 
The  model  is  run  with  a  24  hour  hindcast  and  a  48  hour  forecast.  Personnel  at  the  NAVO 
removed  data  every  6  hours  from  the  first  day  of  the  run  and  stitched  these  Tau’s  together 
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into  a  six  month  long  dataset.  The  area  is  from  23.5°N  to  30.5°  N  and  47.5°E  to  57.5°E. 
The  horizontal  resolution  of  the  operational  model  was  27  km.  The  data  was  retrieved  in 
six  hour  intervals  from  Eebruary  through  July  2006.  Each  six  hour  field  contained  both  U 
and  V  components  of  10  m  winds,  commonly  referred  to  as  surface  winds  and  hereafter 
as  just  winds.  There  were  a  few  days  that  this  could  not  be  accomplished  or  the  data 
collected  was  bad.  This  data  was  removed  during  the  data  analysis  process  and  replaced 
with  the  mean  calculated  from  the  rest  of  the  dataset. 

The  maximum  wind  speed  observed  during  the  entire  period  near  the  oil  terminals 
was  10  m/s.  Based  on  climatology,  and  Eigure  7,  this  value  appears  to  be  low.  During 
Shamal  events  winds  can  reach  up  to  50  kts,  or  25  m/s.  Since  several  Shamals  are  likely 
to  have  occurred  over  a  six  month  period  it  is  extremely  unlikely  that  these  events  never 
created  winds  in  excess  of  10  m/s.  The  fact  that  the  models  do  not  include  wind  gusts  is 
likely  a  factor.  However,  this  would  only  remove  the  more  extreme  events  and  winds  of 
greater  than  10  m/s  are  likely  to  have  occurred.  The  reason  values  greater  than  10  m/s  are 
still  not  observed  can  be  explained  by  the  temporal  resolution.  The  wind  fields  are 
extracted  every  six  hours,  so  events  that  occur  between  these  periods  are  not  observed  in 
model  outputs.  This  limiting  factor  must  be  taken  into  account  during  operational 
planning  and  the  examination  of  the  model  data  in  this  thesis. 

B,  DATA  PROCESSING 

Once  the  data  had  been  collected  and  formatted  for  transfer,  it  was  sent  to  the 
Naval  Postgraduate  School  for  analysis  and  the  currents  were  also  sent  to  ASA.  Once 
ASA  received  the  currents  data,  they  formatted  them  to  be  compatible  with  ChemMap^'^ 
and  sent  the  formatted  data  on  to  Naval  Postgraduate  School  for  use  in  the  thesis.  Once 
the  datasets  were  received,  analysis  and  processing  began.  The  primary  method  of 
analyzing  the  data  was  using  different  methods  available  within  MATrix  EABoratory 
(MATLAB).  MATEAB  is  a  mathematical  software  package  developed  by  Mathworks. 

The  mathematical  and  statistical  analysis  explained  earlier  will  be  made  in  this 
section.  Of  note,  EOE  calculations  become  increasingly  computationally  expensive  as  the 
dimensions  of  the  variables  get  larger.  Beyond  very  small  matrixes  computers  are 
required  to  calculate  EOEs.  As  with  the  data  in  this  study,  the  computer  sometimes 
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becomes  the  limiting  factor.  Subsequently,  sometimes  the  dimensions  of  the  dataset, 
either  spatially,  temporally  or  both,  had  to  be  decreased  in  order  to  complete  the 
calculation.  Since  the  data  provided  for  currents  was  the  larger  and  therefore  more 
difficult  dataset  to  analyze,  most  of  this  section  will  focus  on  currents.  The  wind  dataset 
was  analyzed  with  few  differences  from  that  of  currents.  The  differences  between  winds 
and  currents  will  be  discussed  after  all  methods  for  analyzing  currents  have  been 
mentioned. 

Since  the  current  fields  included  both  data  in  the  water  and  place  holders  for  land, 
the  first  step  in  the  analysis  was  to  determine  which  data  was  associated  with  land  and 
which  data  was  associated  with  water,  and  therefore  currents.  Once  this  was 
accomplished  the  data  associated  with  land  could  be  converted  to  a  variable  that  would 
not  interfere  with  calculations.  This  was  done  via  two  methods.  The  first  method  was  to 
assign  the  land  data  a  value  of  ‘NaN’,  short  for  ‘Not  a  Number’.  By  doing  this  step  when 
a  spatial  plot  was  required  no  values  would  be  placed  on  land.  The  other  method  utilized 
was  to  create  an  index  matrix  of  all  data  associated  with  currents.  By  doing  this,  the  data 
associated  with  currents  could  be  extracted  and  analyzed  without  land  values  interfering. 

After  the  land  data  was  removed,  the  monthly  data  sets  were  combined  to  form 
one  time  series.  Due  to  the  amount  and  size  of  data  collected,  the  dimensions  of  the  data 
set  had  to  be  reduced.  The  first  dimension  that  was  reduced  was  the  overall  spatial 
dimension.  Since  the  area  that  was  being  analyzed  was  determined  to  be  a  low  energy 
area  from  Reynolds  (1993),  the  dataset  dimensions  were  decreased  to  a  box  with 
dimensions: 27.5°  A  to  30.5° N  ;47.5°£’  to  5\° E  .  This  box  was  entirely  in  the  low  energy 
area  thereby  allowing  for  the  study  of  this  regime. 

Focusing  solely  on  the  Northern  Persian  Gulf  greatly  reduced  the  size  of  the  data; 
however  to  complete  the  EOF  calculations  the  dataset  had  to  be  reduced  further.  Instead 
of  analyzing  every  data  point  spatially,  it  was  determined  that  retrieving  every  other  data 
point  still  provided  much  greater  resolution  than  the  wind  data  that  it  would  eventually  be 
compared  to,  yet  still  was  able  to  be  processed  by  the  available  computers.  This  spacing 
remained  fine  enough  to  detect  most  eddies  or  other  phenomena  that  might  occur.  The 
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dataset  was  also  reduced  temporally.  Having  hourly  data  helps  when  studying  tides, 
however  reducing  the  dataset  to  every  third  hour  was  necessary  to  be  able  to  compute  the 
EOFs.  This  reduction  in  the  dataset  was  primarily  used  for  the  EOF  calculations.  For 
example,  when  the  time  series  calculations  were  done  on  a  location  near  Iraq’s  Oil 
Terminals,  the  time  series  included  all  of  the  hourly  datapoints. 

The  next  step  was  to  see  if  there  was  any  data  that  was  obviously  corrupted.  This 
was  accomplished  while  computations  of  EOFs  were  being  run.  While  analyzing  the 
EOFs  there  was  a  period,  illustrated  in  Figure  12,  which  shows  extremely  different  than 
the  rest  of  the  data.  This  data,  corresponding  to  1 8  May,  was  determined  to  be  erroneous 
and  subsequently  replaced  by  the  mean  of  the  remaining  data.  The  mean  was  used  due  to 
the  fact  that  in  the  EOF  calculations  the  mean  would  be  removed  and  therefore  the  data 
would  essentially  be  zero.  It  was  also  noted  by  NAVO  scientists  that  the  data  associated 
with  23  June  was  missing  due  to  computer  problems  on  that  day.  The  mean  of  the  data 
was  inserted  into  this  date  for  the  reasons  discussed  earlier.  These  two  periods  were 
removed  from  consideration  when  examining  times  to  run  the  chemical  models. 


Figure  12.  Time  Series  of  EOF  Coefficients  for  the  currents  of  the  Northern  Persian 

Gulf  on  a  semi-log  y  scale. 
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Once  the  dataset  had  been  converted  into  a  form  that  was  able  to  be  examined  and 
all  days  were  filled  with  either  aetual  data  or  replaced  with  the  mean,  analysis  of  the  data 
eould  begin.  Several  goals  associated  with  the  analysis  ineludes:  (1)  examine  a  time 
series  of  a  point  near  the  oil  terminals  to  extract  times  for  the  model  runs,  (2)  determine 
how  many  EOF  modes  were  required  to  eapture  the  majority  of  the  Northern  Persian 
Gulfs  eharaeteristics,  and  (3)  examine  the  data  in  a  spatial  environment  to  get  a  feel  for 
primary  the  amplitudes  and  directions  of  currents. 

C.  WINDS 

1.  Monthly  Mean 

Analysis  of  the  wind  began  with  examining  the  mean  wind  over  the  entire  time 
and  entire  Persian  Gulf  and  progressed  to  a  more  detailed  analysis  of  each  month’s 
characteristics.  The  mean  wind  for  the  entire  time.  Figure  13,  showed  a  distinct  flow 
along  the  major  axis.  The  predominantly  northwesterly  flow  resembles  the  Shamal 
(Perrone,  1979).  As  the  flow  approached  the  Straits  of  Hormuz,  it  veered  towards  Iran’s 
coastline.  Flow  in  the  Gulf  of  Oman  also  tends  towards  Iran’s  eoast.  A  likely  eause  for 
the  flow  to  turn  north  in  both  instanees  is  the  al-Hajar  mountain  range  along  the 
Musandam  Peninsula.  This  mountain  range  would  tend  to  obstruct  flow  and  direct  it 
along  the  mountain  range.  The  mean  flow  in  Iran  tends  to  flow  along  the  Zagros 
Mountains.  Mean  flow  to  the  southwest  of  the  Persian  Gulf  is  weak,  with  a  general 
direction  to  the  south.  As  seen  in  Figure  13,  the  mean  wind  speed  is  less  than  5  m/s 
throughout  the  Persian  Gulf  Near  Iraq’s  oil  terminals,  the  mean  wind  speed  is  even 
weaker,  around  2  -3  m/s. 
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Figure  13.  Mean  wind  vectors  of  the  entire  dataset  (Feb  -  Jul). 


When  looking  at  the  months  individually  it  became  clear  that  the  characteristics 
and  magnitudes  changed,  sometimes  drastically,  from  month  to  month.  From  examining 
February’s  mean  wind  field,  it  is  easy  to  see  there  is  a  line  that  starts  near  2S° N  50° E 
and  extends  along  the  major  axis  until  it  reaches  the  shore  (Figure  14a).  This  is  a  line  of 
bifurcation  between  flow  turning  towards  the  Arabian  Peninsula  and  the  Asian  Continent. 
There  is  also  distinct  area  of  weak  winds  surrounding  Iraq’s  oil  terminals.  As  the  winds 
progress  down  the  major  axis  of  the  Persian  Gulf,  their  magnitudes  tend  to  increase.  This 
increase  may  be  due  to  the  decreased  friction  resulting  from  the  winds  traveling  from 
land  to  water  and  enhancement  downslope  flow  off  the  Zagros  Mountains  to  the  east. 
Relatively  same  magnitudes  were  seen  in  March;  however  the  flow  was  slightly  different. 
The  region  of  weak  winds  near  the  oil  terminals  has  increased  in  size,  and  the  primary 
direction  of  those  winds  are  now  across  the  major  axis  (Figure  14b).  The  line  of 
bifurcation  has  also  migrated.  It  now  extends  from  the  northern  tip  of  Qatar  along  the 
major  axis  to  the  shore.  Flow  north  of  Qatar  is  predominantly  along  the  major  axis  with 
only  small  deviations  away  from  this  continuity. 
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By  April  the  area  near  the  oil  terminals  has  inereased  further  to  the  south  and  an 
anti-oyelonie  pattern  is  elearly  visible  (Figure  14e).  Besides  the  inerease  in  the  area  to  the 
north,  the  winds  over  the  Persian  Gulf  showed  remarkably  little  change  from  March  to 
April.  However,  this  relatively  constant  regime  disappeared  in  May.  As  illustrated  in 
Figure  14  d,  the  magnitudes  throughout  the  region  increased  significantly,  although  the 
pattern  only  changed  in  the  north.  There  was  no  longer  an  area  of  weak  winds 
surrounding  the  oil  terminals.  The  winds  predominantly  flowed  along  the  major  axis  and 
either  over  the  Musandam  Peninsula  or  into  the  Straits  of  Hormuz.  If  the  winds  traveled 
into  the  straits,  they  ended  up  drastically  decreasing  once  the  encountered  lands,  and 
therefore  the  Zagros  Mountains.  Effects  from  the  Zagros  Mountains  are  also  felt 
throughout  the  eastern  side  of  the  Persian  Gulf.  It  is  evident  that  along  this  side  the 
Zagros  Mountains  due  on  of  two  things.  First,  they  could  block  the  flow  and  subsequently 
direct  it  along  the  east  coast  of  the  Persian  Gulf.  The  other  effect  is  flow  coming  off  the 
Zagros  Mountains.  This  flow  converges  with  the  flow  in  the  Persian  Gulf  thus  enhancing 
the  already  stronger  flow. 

By  June,  the  presence  of  the  Southwest  Monsoon  in  the  Arabian  Sea  is  felt  in  the 
Gulf  of  Oman.  Flow  is  no  longer  away  from  the  Musandam  Peninsula,  but  in  from  the 
Arabian  Sea  (Figure  14e).  This  flow  is  not  the  only  significant  change.  In  the  region 
around  the  oil  terminals,  which  previously  was  an  area  of  weak  winds,  there  is  moderate 
flow  down  the  major  axis.  The  line  of  bifurcation  no  longer  exists  south  of  Qatar,  being 
replaced  by  a  new  area  of  weak  winds.  Flow  on  the  eastern  edge  of  the  Persian  Gulf  show 
a  strong  area  of  convergence.  The  coastline  now  marks  where  downslope  flows  from  the 
Zagros  Mountains  oppose  moderate  flow  from  the  northern  Persian  Gulf.  With  nowhere 
to  go,  this  area  of  convergence  is  also  an  area  of  upward  vertical  motion. 

The  majority  of  flow  in  July  is  similar  to  that  seen  in  June.  The  Gulf  of  Oman  still 
feels  the  influence  of  the  Southwest  Monsoon  and  the  flow  in  the  northern  Persia  Gulf  is 
still  moderate  (Figure  14f).  What  differs  is  the  area  of  weak  flow  just  south  of  Qatar  and 
the  direction  of  flow  in  the  north.  The  region  of  weak  flow  has  become  virtually 
nonexistent  with  the  majority  of  flow  now  moving  towards  the  Arabian  Peninsula.  The 
direction  of  flow  to  the  north  in  no  longer  along  the  major  axis,  instead,  it  has  a  distinct 
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across  axis  portion.  This  equates  to  flow  reaehing  the  east  eoast  while  still  in  the  northern 
region.  Thus  the  feteh  is  markedly  redueed  and  there  is  less  time  for  the  magnitude  to 
build. 
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Mean  Wind  Vectors  for  FEB 


Mean  Wind  Vectors  for  MAR 
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b) 


Mean  Wind  Vectors  for  APR 


Mean  Wind  Vectors  for  MAY 
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c) 


d) 


Mean  Wind  Vectors  for  JUN 


Mean  Wind  Vectors  for  JUL 
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e) 


Figure  14. 


f) 

Mean  monthly  wind  vectors. 
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2,  Anomalies 

Anomalies  were  primarily  examined  over  the  entire  record  length  and  in  the  form 
of  analyzing  the  first  three  EOF  modes.  As  mentioned  earlier,  to  get  a  better  picture  of  the 
low  energy  environment  examined  in  this  thesis,  EOFs  were  conducted  only  in  the 
northern  Persian  Gulf  If  the  variance  calculated  during  the  EOF  process  is  summed  up 
and  then  each  variance  is  divided  by  that  sum,  a  percent  of  total  variance  is  obtained. 

var 

%var= - ^  (7.1) 

1=1 

Since  the  variance  was  obtained  from  computer  based  EOF  calculations,  there 
were  a  finite  number  of  modes  that  could  be  derived  from  the  data.  The  wind  data 
allowed  for  a  total  of  300.  From  these  calculations  it  was  determined,  as  illustrated  in 
Table  1  below,  that  the  first  3  modes  encompassed  88  %  of  the  variance.  By  the  fourth 
mode  the  percent  of  the  total  variance  has  decreased  drastically  and  had  less  significance. 


EOF  mode 

Percent  of  Total  Variance 

Cumulative  Percent 

First  Mode 

61.1860 

61.1860 

Second  Mode 

20.1800 

81.3660 

Third  Mode 

6.6906 

88.0566 

Fourth  Mode 

3.2760 

91.3326 

Fifth  Mode 

1.9524 

93.2850 

Table  1.  Percent  of  Variance  for  the  first  5  EOF  modes  of  Wind  data 

Based  on  the  percent  variance,  the  first  three  EOF  modes  were  analyzed  to 

determine  likely  forcing.  Although  these  modes  were  examined,  the  wind  field  depicted 

in  them  is  significantly  less  than  what  was  examined  in  the  mean  wind  vector  discussion. 

This  is  due  to  the  fact  that  EOFs  looked  at  the  perturbations.  To  get  an  idea  of  how  large 

of  winds  the  EOFs  produce,  take  the  representative  vector  from  Figure  15  and  multiply  it 

by  the  magnitude  scale  in  Figure  16.  When  this  is  done  for  EOF  Modes  1,  2,  and  3,  the 

scales  are  shown  to  be  around  2  m/s,  0.8  m/s,  and  0.8  m/s,  respectively. 
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EOF-1  has  a  relatively  strong  uniform  eastward  flow  (Figure  15).  The  angle  that 
this  EOF  is  multiplied  by  is  predominantly  close  to  either  0  degrees  or  close  to  180 
degrees  (Figure  16).  This  indicates  the  primary  flow  other  than  the  mean  is  across  the 
gulf  When  the  angle  is  close  to  0  degrees  the  flow  is  like  in  Figure  15,  for  180  degrees 
the  flow  is  in  the  reverse  direction.  This  flow  likely  corresponds  to  passage  of  fronts  that 
would  allow  for  higher  winds  across  the  gulf  The  winds  flowing  down  the  gulf  are 
already  accounted  for  in  the  mean.  This  conclusion  makes  more  sense  when  the 
magnitudes  of  the  winds  are  also  examined.  During  the  winter,  February  into  March, 
there  are  higher  magnitudes,  occurring  every  few  days.  The  magnitudes  then  decrease 
until  the  next  event.  This  timeframe  is  consistent  with  the  findings  of  Perrone,  1979. 

The  next  mode,  EOF-2,  has  a  cyclone  centered  near30°A^  48.75°F'.  A  cull  is 
located  at  the  edge  of  the  area,  27°  A  50.5°F'  (Figure  17).  The  Magnitude  and  direction  of 
the  coefficients  are  harder  to  utilize  to  determine  the  rationale  behind  the  circulation.  This 
is  because  the  magnitude  does  not  vary  in  set  intervals  and  does  not  vary  between  only 
two  main  directions  (Figure  18).  The  magnitudes  have  no  set  interval  and  cycle  between 
larger  and  smaller  magnitude  repeatedly  followed  by  period  of  small  magnitudes.  The 
direction  tends  to  continually  change.  Based  on  both  the  previously  mentioned  features 
and  the  fact  that  the  pattern  in  cyclonic  in  nature,  this  likely  the  period  intensification  of  a 
thermal  low  just  landward  of  the  oil  terminals.  Thermal  Fows  are  known  events 
throughout  the  region,  with  the  strongest  Thermal  Fows  found  in  Saudi  Arabia  and  in 
Iraq  (Ackerman  and  Cox,  1982).  The  Thermal  Fows  normally  occur  onshore,  but  near  the 
water.  This  is  due  to  the  increased  pressure  difference  created  by  the  heating  of  the  land 
while  the  water  remains  relatively  constant. 

In  EOF-3  there  is  a  line  of  bifurcation  across  the  Persian  Gulf,  splitting  flow 
curving  to  the  north  and  flow  curving  to  the  south  (Figure  19).  When  looking  at  the 
amplitudes  and  directions  for  EOF-3,  Figure  20,  the  primary  directions  are  near  0  degrees 
and  near  270  degrees.  Magnitudes  are  relatively  small  with  the  times  of  large  magnitudes 
correspond  to  the  same  periods  as  those  of  EOF-1.  These  times  also  correspond  to 
migratory  lows.  Since  EOF-3’s  magnitudes  are  close  to  half  that  of  EOF-1,  EOF-3  would 

correspond  to  residual  influences  of  transiting  low  pressure  systems. 
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Figure  15.  EOF-1  for  entire  wind  dataset  in  the  Northern  Persian  Gulf. 
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Figure  16.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  wind  dataset  in  the 

Northern  Persian  Gulf 


(Theta  Mode  1) 


SWAFS  Magnitude/Theta  Time  Series:  (Magnitude  Model) 
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EOF  of  Wind  (2nd  Mode) 
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Figure  17.  EOF-2  for  entire  wind  dataset  in  the  Northern  Persian  Gulf. 
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Figure  18.  EOF-2  amplitude  (top)  and  direction  (bottom)  for  entire  wind  dataset  in  the 

Northern  Persian  Gulf 
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EOF  of  Wind  (3rd  Mode) 


SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  Mode3) 


Figure  20.  EOF-3  amplitude  (top)  and  direction  (bottom)  for  entire  wind  dataset  in  the 

Northern  Persian  Gulf 
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3,  Spatial  Examining  Raw  Data 

While  the  mean  is  widely  used,  when  examining  a  possible  oscillating  pattern  to 
determine  appropriate  times  to  run  the  models,  the  mean  gives  only  a  partial  view. 
Maxima  at  either  end  of  the  spectrum  are  not  shown.  Therefore,  to  visually  examine  the 
data  and  focus  the  research  on  important  periods  the  mean  was  not  used.  To  accomplish 
this,  looking  at  the  raw  data  in  time  series  format  was  beneficial.  This  was  done  in  two 
steps.  The  first  step  also  accomplished  goal  two.  Specifically,  the  U  and  V  components 
were  converted  to  polar  coordinates.  This  allowed  the  magnitudes  and  directions  to  be 
analyzed  separately.  The  magnitudes  were  a  significant  factor  when  determining  model 
run  times.  Some  of  the  criteria  for  the  scenarios  included  having  a;  (1)  minimal  wind  and 
weak  current  event,  (2)  maximum  wind  and  weak  current  event,  (3)  maximum  wind  and 
maximum  current  event  and  (4)  looking  for  some  period  that  a  weather  phenomena 
occurred.  It  could  easily  be  seen  when  the  wind  criteria  was  met  by  looking  at  a  time 
series  of  the  magnitude  of  winds  (Figure  21). 


to  maximum  value  for  the  definition  of  weak  winds,  while  the  blue  line  corresponds  to 
the  minimum  value  for  strong  winds,  in  this  thesis. 
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Although  when  the  term  maximum  wind  event  are  used  gale  force  or  greater 
winds  are  thought  of,  in  the  low  energy  region  of  the  northern  Persian  Gulf  this  is  not  the 
case.  An  extreme  high  wind  event  was  defined,  for  the  purpose  of  this  study,  as  any  time 
period  that  consistently  exceed  5  m/s.  On  the  opposite  spectrum,  a  low  wind  event  was 
defined  as  any  time  period  that  consistently  was  below  3.5  m/s.  Maximum  wind  events 
occurred  primarily  later  in  the  time  period,  with  the  strongest  event  occurring  in  early 
July.  Weak  wind  events  happened  less  frequently,  with  the  primary  event  occurring  in 
late  February. 

D.  CURRENTS  FOR  THE  NORTHERN  GULF 

The  same  concept  utilized  for  the  wind  analysis  was  applied  to  the  low  energy 
currents  regime  of  interest,  the  northern  Persian  Gulf.  Although  the  results  from  mean 
currents  calculations  were  the  same  when  calculating  the  entire  gulf  or  just  the  northern 
Persian  Gulf,  examining  the  northern  region  allowed  for  a  better  spatial  resolution  to  be 
displayed. 

The  mean  vector  over  the  entire  time  period  and  the  monthly  mean  vector  plots 
appear  to  show  a  turbulent  flow  throughout  the  northern  region  (Figures  22  and  23). 
Keeping  in  mind  that  the  flow  in  this  region  is  primarily  tidally  driven,  it  could  be 
expected  that  the  mean  should  be  zero,  or  no  flow.  Indeed  the  mean  flow  is  actually  very 
small  in  magnitude.  Notice  the  reference  arrow  is  10  cm/s  and  most  of  the  flows  are  less 
than  the  reference.  The  mean  flow  shows  preferred  paths  of  flow  and  not  what  happens  at 
any  one  time  throughout  the  given  period. 

The  mean  currents  for  the  entire  time.  Figure  22,  show  flow  north  along  the 
Iranian  coast,  across  the  gulf  and  down  along  the  Arabian  Peninsula  and  central  gulf  This 
pattern  extends  up  to  29.5°  N  .  North  of  29.5°  N  was  a  region  of  primarily  slow  southward 
flow. 
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Vector  Mean  of  Currents 


A8°£  49°E  50“E  51°E 

Longitude 

Figure  22.  Mean  current  vectors  of  the  entire  dataset  (Feb  -  Jul)  Arrow  shows  a  10 

cm/s  current. 


When  looking  at  the  months  individually  it  became  clear  that  the  characteristics 
and  magnitudes  changed  drastically  from  the  first  few  months  to  the  last  months 
collected.  However,  every  month  contained  a  current  running  southward  near  the  major 
axis  of  the  gulf. 

From  examining  February’s  mean  current  field,  the  flow  throughout  the  region  is 
weak  and  less  organized  then  the  mean  for  the  entire  time  span  (Figure  23  a).  There  is  a 
weak  flow  up  the  Arabian  Peninsula.  This  flow  then  turns  towards  the  central  region  and 
combines  with  flow  coming  from  the  north  and  east.  This  combined  flow  forms  the 
strongest  current  for  February.  This  combined  flow  is  the  characteristic  southward  flow 
observed  in  every  month  down  the  major  axis.  Flow  north  of  29°Atends  to  be  very 
weak,  less  than  2.5  cm/s. 

Approximately  the  same  magnitude  was  seen  in  March;  while  the  character  of  the 

flow  has  changed  some  (Figure  23b).  The  characteristic  southward  flow  observed  in 
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every  month  down  the  major  axis  begins  a  little  further  north.  The  flow  up  the  Arabian 
Peninsula  has  dissipated  being  replaeed  by  weak  southward  flow  originating  in  the 
extreme  northwest.  There  is  a  commencement  of  flow  up  the  Iranian  coast  seen  in  the 
mean  flow  for  the  entire  period;  although  it  is  not  yet  well  organized  Flow  north  of  29°  N 
remains  weak,  but  is  becoming  more  organized.  Flow  across  the  gulf  is  commencing  in 
the  northern  region.  The  flow  near  the  oil  terminals  remains  weak  and  southward.  The 
mean  flow  in  April,  Figure  23c,  sees  a  continuation  of  the  same  trends  seen  in  March. 
Flows  along  the  Arabian  Peninsula  and  up  the  Iranian  coast  have  intensified.  The 
intensification  and  general  flow  of  April  makes  it  the  most  similar  to  the  mean  flow.  The 
main  difference  between  April  and  the  entire  mean  flow  is  the  central  region.  In  April  this 
region  is  still  unorganized  and  several  eddy-like  formations  are  observed;  whereas  in  the 
mean  flow  for  the  entire  period,  the  central  region  is  well  defined  and  flow  primarily 
along  the  major  axis  to  the  south. 

By  May,  the  overall  magnitudes  of  the  mean  flow  has  increased  significantly 
(Figure  23  d).  Although  the  characteristic  southward  flow  is  still  observed  in  May,  it  is 
partially  hidden  by  the  fact  that  most  of  the  flow  south  of  29.5° N  is  in  the  same  general 
direction  and  has  the  same  magnitude.  It  is  evident  there  is  still  flow  up  the  Iranian  coast; 
however,  this  flow  is  now  masked  by  the  predominant  southerly  and  across  gulf  flows. 

The  mean  flow  in  June  is  similar  to  the  entire  periods  mean  flow  (Figure  23  e). 
The  major  difference  is  in  the  relative  strength  of  the  flow.  There  is  a  predominant  flow 
up  the  Iranian  coast,  across  the  gulf,  and  back  down  along  the  Arabian  Peninsula. 
Although  the  magnitude  for  the  entire  time  period  tends  to  be  less  than  10  cm/s,  the  mean 
for  June  tends  to  be  at  or  above  that.  Flow  in  the  extreme  northern  region,  near  the  oil 
terminals,  remains  low.  When  examining  the  mean  flow  for  July,  Figure  23f,  a 
predominant  flow  to  the  south  is  observed.  A  few  eddy-like  features  near  the  Iranian 
coast  are  also  observed,  but  the  primary  change  is  in  the  extreme  north.  It  is  here  that  the 
flow  has  started  to  increase.  No  longer  is  there  just  a  weak  flow,  the  currents  are  now  to 
the  south  at  around  5  cm/s.  This  is  a  major  change  from  the  near  no  flow  in  the  previous 
months. 
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c) 


d) 


e) 


f) 

Figure  23.  Mean  monthly  eurrent  vectors. 
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Anomalies  were  examined  over  the  entire  record  length  and  in  the  form  of 
analyzing  the  first  three  EOF  modes.  As  mentioned  earlier,  to  get  a  better  picture  of  the 
low  energy  environment  examined  in  this  thesis,  EOFs  were  conducted  only  in  the 
northern  Persian  Gulf. 

The  same  method  used  to  calculate  percent  variance  for  winds,  as  illustrated  in 
equation  (7.1),  was  applied  to  currents.  Since  the  variance  was  obtained  based  EOF 
calculations,  there  were  a  finite  number  of  modes  that  could  be  derived  from  the  data. 
The  currents  data  allowed  for  a  total  of  over  500  modes.  The  reason  for  more  modes  in 
currents  than  winds  is  due  to  the  higher  spatial  and  temporal  resolution  of  the  currents. 
The  modes  for  currents  were  truncated  at  500  for  several  reasons.  The  primary  reason 
was  at  mode500  the  percent  of  the  mode  one  variance  was  less  than  one  ten  thousandth  of 
a  percent,  and  therefore  all  values  thereafter  would  be  negligible.  Another  less  significant 
reason  was  processing  time.  From  these  calculations  it  was  determined,  as  illustrated  in 
Table  2  below,  that  the  first  3  modes  encompassed  91  %  of  the  variance.  By  the  fourth 
mode  the  percent  of  the  total  variance  has  decreased  drastically  and  had  less  significance. 


EOF  mode 

Percent  of  Total  Variance 

Cumulative  Percent 

First  Mode 

69.3650 

69.3650 

Second  Mode 

16.4190 

85.7840 

Third  Mode 

5.3307 

91.1147 

Fourth  Mode 

1.6572 

92.7719 

Fifth  Mode 

1.1662 

93.9381 

Table  2.  Percent  of  Variance  for  the  first  5  EOF  modes  of  Current  data 


Based  on  the  percent  variance,  the  first  three  EOF  modes  were  analyzed  to 
determine  likely  forcing.  Although  these  modes  were  examined,  the  current  field  depicted 
in  them  is  significantly  less  than  what  was  examined  in  the  mean  wind  vector  discussion. 
This  is  due  to  the  fact  that  EOFs  looked  at  the  perturbations.  To  get  an  idea  of  how  large 
of  currents  the  EOFs  produce,  take  the  representative  vector  from  Figure  24  and  multiply 
it  by  the  magnitude  scale  in  Figure  25.  When  this  is  done  for  EOF  Modes  1,  2,  and  3,  the 
scales  are  shown  to  be  around  100  cm/s,  50  cm/s,  and  25  cm/s,  respectively. 
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After  examining  the  magnitudes  and  direetions  of  EOF-1  (Figure  25),  EOF-2 
(Figure  27)  and  EOF-3  (Figure  29),  it  is  apparent  that  tidal  influenees  are  evident 
throughout  the  first  three  modes.  Although  the  magnitudes  were  of  different  seales,  eaeh 
have  peaks  and  lulls  at  approximately  the  same  time  periods.  On  a  shorter  seale,  the 
maximums  are  always  immediately  followed  by  a  minimum.  The  direetion  of  the  first 
two  modes  changed  primarily  by  180  degrees.  The  two  distinct  periods  found  in  the 
magnitudes  of  the  first  three  EOFs  and  the  directional  shift  in  the  first  two  EOFs  indicates 
there  is  a  strong  tidal  influence  in  these  modes.  Further  examination  was  done  to  try  and 
determine  the  other  primary  influence  within  the  EOFs.  If  no  other  major  stimulus  is 
observed,  further  examination  of  the  tidal  influences  will  be  discussed.  While  the  tidal 
influences  were  readily  apparent  in  all  three  modes,  there  is  a  possible  amphidrome 
located  within  the  area  of  low  energy.  Figure  8,  which  would  make  it  impossible  to  filter 
out  the  tidal  influences  from  the  EOFs.  An  attempt  to  isolate  the  tidal  influences  was 
made  using  a  low  pass  fdter;  however,  the  influences  were  still  apparent. 

EOF-1  has  predominantly  southwesterly  flow  (Figure  24).  After  comparing  the 
magnitudes  to  the  mean  wind  magnitudes,  it  was  determined  that  the  flow  corresponding 
to  EOF-1  has  a  strong  wind-driven  component.  The  direction  of  flow.  Figure  25,  was 
predominantly  near  0  degrees  and  near  180  degrees.  When  the  angle  was  near  180 
degrees,  it  would  slowly  decrease  from  where  it  started  for  a  period  of  time,  normally 
decreasing  approximately  45  degrees.  When  comparing  the  magnitudes  and  spatial 
picture  of  EOF-1  from  currents.  Figures  24  and  25,  to  the  wind  magnitude  time  series  and 
mean  wind  vector  fields.  Figures  21  and  22,  as  direction  correlation  can  easily  be 
established.  When  wind  magnitudes  increase  so  do  the  EOF-1  magnitudes.  Although  the 
direction  is  deceiving,  the  fact  that  EOF-1  tend  to  be  at  180  degrees  when  the 
corresponding  magnitudes  are  maximum,  reverses  the  direction  of  interest  and 
corresponds  to  wind-driven  currents.  The  gradual  shift  of  the  direction  in  EOF-1 
corresponds  to  the  time  lag  in  the  direction  due  to  wind  forcing 

The  mechanism  behind  EOF-2  is  much  more  apparent  (Figures  26  and  27).  The 
period  from  one  angle  back  to  the  same  angle  is  around  14  days.  It  actually  varies 
between  13  and  16  days.  The  direction  cycles  between  just  under  90  degrees  to  just  under 
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270  degrees.  The  amplitude  has  a  cyclic  maximum  on  the  same  time  scale  as  the  angle, 
where  a  maximum  corresponds  to  a  corresponding  large  angle.  When  looking  at  the 
vector  plot  of  EOF-2,  there  is  a  distinct  cyclonic  flow.  All  of  the  previously  mentioned 
factors  point  towards  EOF-2  being  tidally  driven.  As  mentioned  previously  all  the  modes 
had  a  tidal  influence;  however,  EOF-2  appeared  to  have  no  other  dominating  influences. 

EOF-3  shows  influences  from  several  different  factors  (Figures  28  and  29).  The 
magnitudes  cycle  on  the  same  time  scale  as  spring  and  neap  tides.  The  Direction, 
however,  does  not  illustrate  a  dominant  tidal  influence.  In  fact  there  is  no  dominant  cycle 
established  within  the  direction  time  series.  The  vector  plot  of  EOF-3  also  does  not  have 
a  strong  tidal  influence.  However,  there  are  two  main  areas  in  this  EOF,  one  on  each  side 
of  the  gulf.  The  eastern  side  has  a  line  of  convergence  running  along  the  Iranian  coast; 
while  the  western  side  has  a  parallel  line  of  divergence.  The  cause  of  the  convergence  and 
divergence  is  unknown.  However,  with  the  percent  variance  of  EOF-3  equal  to  around 
5%  and  the  magnitude  corresponding  to  a  maximum  of  25  cm/s  this  mode  has  little 
significance  compared  to  modes  1  and  2.  The  fact  that  there  is  not  a  single  dominating 
factor  and  the  factors  that  can  be  distinguished  are  accounted  for  in  the  previous  EOFs 
and  the  mean  current,  further  diminish  the  relevance  of  this  mode.  In  the  end,  85%  of  the 
variance  was  accounted  for  in  the  first  two  modes  and  some  of  the  influences  were 
established  for  EOF-3,  but  determining  all  of  the  possible  influences  for  EOF-3  was  left 
for  future  consideration. 
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Longitude 

Figure  24.  EOF-1  for  entire  eurrents  dataset  in  the  Northern  Persian  Gulf 


SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  Model) 


Figure  25.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  currents  dataset  in 

the  Northern  Persian  Gulf. 
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Magnitude 


Longitude 

Figure  26.  EOF-2  for  entire  currents  dataset  in  the  Northern  Persian  Gulf 


SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  Model) 


(Theta  Mode  2) 


the  Northern  Persian  Gulf. 
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Figure  28.  EOF-3  for  entire  currents  dataset  in  the  Northern  Persian  Gulf 


SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  ModeS) 

n  ^  ^  ^  1^ 


(Theta  Mode  3) 


O  ^  ^  ..gT  V  <1,^  .01^  ^  ^  ^  ^  ^  <0  ^  ^  ^ 

Figure  29.  EOF-3  amplitude  (top)  and  direction  (bottom)  for  entire  currents  dataset  in 


the  Northern  Persian  Gulf. 
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3,  Spatial  Examining  Raw  Data 

As  discussed  in  the  winds  seetion,  when  examining  a  possible  oseillating  pattern 
to  determine  appropriate  times  to  run  the  models,  the  mean  gives  only  a  partial  view.  The 
same  logie  utilized  in  examining  winds  was  also  used  to  examine  eurrents.  To  further 
examine  the  eurrents,  looking  at  the  raw  data  in  time  series  format  was  benefieial.  This 
was  done  in  two  steps.  The  first  step  also  aeeomplished  goal  two.  Speeifieally,  the  U  and 
V  eomponents  were  eonverted  to  polar  eoordinates.  This  allowed  the  magnitudes  and 
direetions  to  be  analyzed  separately.  The  magnitudes  were  a  signifieant  faetor  when 
determining  model  run  times.  As  diseussed  in  the  winds  seetion,  some  of  the  eriteria  for 
the  seenarios  ineluded  having  a:  (1)  minimal  wind  and  weak  eurrent  event,  (2)  maximum 
wind  and  weak  eurrent  event,  (3)  maximum  wind  and  maximum  eurrent  event  and  (4) 
looking  for  some  period  that  a  weather  phenomena  oeeurred.  Most  of  these  eriteria  eould 
easily  be  seen  when  looking  at  a  time  series  of  the  magnitude  of  eurrents  (Figure  30)  and 
eomparing  it  to  that  of  winds  (Figure  21). 

Although  when  the  term  maximum  eurrent  event  is  used  several  m/s  eurrents 
often  are  thought  of,  in  the  low  energy  region  of  the  northern  Persian  Gulf  this  is  not  the 
ease.  A  high  eurrent  event  was  defined  as  any  event  where  eurrents  repeatedly  exeeeded 
40  em/s,  with  a  low  eurrent  event  being  one  that  remained  below  25  ern/s.  The  term 
repeatedly  is  used  when  defining  eurrents  due  to  the  eyelie  nature  of  eurrents.  Several 
periods  of  maximum  eurrents  were  found,  primarily  in  late  May  and  through  July. 
Periods  of  weak  eurrents  were  found  through  out  the  time  period,  ineluding;  (1)  the  last 
week  in  February,  (2)  the  first  week  in  April,  and  (3)  the  first  week  in  July. 
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Timeseries  of  the  Amplitude  of  Currents 


Figure  30.  Time  Series  of  the  Magnitude  of  Currents  (in  em/s).  The  red  line 
eorresponds  to  maximum  value  for  the  definition  of  low  eurrents,  while  the  blue  line 
eorresponds  to  the  minimum  value  for  high  eurrents,  in  this  thesis. 


E,  EXTRACTING  TIMES  FOR  CHEMICAL  MODEL  RUNS 

Comparing  the  wind  and  eurrents  datasets  was  a  relatively  straightforward 
proeess.  This  was  aeeomplished  by  ereating  a  MATLAB  file  that  overlaid  the  veetor 
wind  field  over  a  eurrents  field  on  a  temporal  domain  of  every  six  hours.  An  example  of 
this  is  given  in  Figure  31.  Onee  this  was  aeeomplished,  a  slideshow  was  ereated  and  run. 
This  allowed  a  visual  inspeetion  of  the  wind  and  eurrents  fields.  Dates  were  then 
extracted  that  met  model  scenario  criteria.  These  dates  were  then  compared  to  the  time 
series  created  earlier  of  winds  and  of  currents  (Figure  21  and  30).  Once  the  chosen  times 
were  confirmed  in  both  the  time  series  analysis  and  temporal  slideshow,  model  runs  were 
commenced. 
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Vector  Comparison:  nliids  &  current(ni/s)  in  SOU  2006,  14  May,  12  CTC,  16  LST 


51  °E 


49°E 
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Figure  31.  Comparison  of  Wind  and  Currents  fields  from  14  May  at  12Z. 
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VIII.  CONNECTION  BETWEEN  PHYSICAL  AND  CHEMICAL 

MODELS 


The  process  of  taking  the  output  of  the  physical  models  and  using  it  as  an 
environmental  input  for  the  chemical  models  was  long  and  arduous.  This  section  will 
detail  the  process  by  which  COAMPS  output  was  formatted  to  be  ingested  into  HPAC.  It 
will  also  describe  how  SWAPS  data  was  formatted  and  ingested  into  CHEMMAP™. 
Lastly,  it  will  describe  the  process  by  which  some  of  the  COAMPS  data  was  ingested  into 
CHEMMAP™. 

To  ingest  COAMPS  wind  fields  into  HPAC  the  fields  must  be  converted  to  a 
format  HPAC  recognizes.  This  process  was  done  in  Matlab.  To  start  the  process  a  header 
that  followed  the  format  in  Table  3  was  created.  The  lines  following  this  header  included 
the  data  from  COAMPS  for  the  times  needed.  Several  files  were  made,  one  for  each 
model  runs  time,  using  this  technique.  These  files  were  then  opened  using  the  weather 
file  editor  and  ingested  the  new  file  back  into  HPAC.  The  wind  was  then  utilized  as  a 
surface  wind  over  the  region. 


n 

CREATOR: 

Universal Weather Reader 

(Chuckles/ 172 . 

.20.69. 105) 

n 

DATE  : 

2006 

-10 

-03 

21:02:35  UTC 

# 

EDITED : 

yes 

n 

REFERENCE : 

agl 

u 

ANALYSIS : 

2006 

02 

01 

00  .  OO 

n 

START: 

2006 

02 

oi 

OO  .  OO 

n 

END  : 

2006 

02 

02 

12  .  OO 

u 

TIME REFERENCE ; 

;  utc 

n 

MODE  : 

otos 

all 

SURFACE 

8 

ID  YYYYHHDDHOUR 

LAT 

LON 

Z 

WDIR 

WSPD 

HOURS 

N 

E 

H 

DEG 

H/S 

-9999 

Table  3.  Example  of  HPAC  file  header  required  to  ingest  data  into  HPAC  from 

archived  COAMPS  fields. 


Ingesting  the  data  into  CHEMMAP  was  an  entirely  different  problem.  After 
trying  to  accomplish  it  by  following  the  same  methods  to  ingest  COAMPS  into  HPAC,  it 
was  discovered  that  it  was  impossible  to  do  from  the  Naval  Postgraduate  School.  ASA 
was  very  helpful  and  offered  to  set-up  a  file  transfer  site  and  configure  the  datasets  to  a 
usable  format  for  CHEMMAP™.  This  assistance  was  the  only  way  to  get  archived  data 
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into  the  proper  format  in  a  timely  manner.  However,  if  CHEMMAP™  had  been  run  real¬ 
time,  the  data  eould  have  been  retrieved  from  the  operational  models  directly,  thereby 
eliminating  the  configuration  problem.  Due  to  time  constraints,  the  winds  entire  wind 
field  was  not  able  to  be  converted  to  a  spatially  and  temporally  varying  file.  However,  a 
time  series  of  a  location  near  Iraq’s  oil  terminals  was  able  to  be  processed  and  ingested 
into  CHEMMAP™.  Once  this  field  was  ingested  into  CHEMMAP™  it  was  treated  as  a 
uniform  wind  field,  thereby  providing  a  wind  vector  for  the  entire  region. 

A  drawback  to  the  fields  that  were  obtained  via  the  physical  models  was  in  the 
limited  variables  archived.  Although  the  chemical  models  would  take  more  detailed 
inputs,  the  fields  utilized  were  at  only  one  depth,  or  height.  There  were  several  other 
possible  inputs  into  the  chemical  models  that  could  have  been  utilized  if  the  fields  were 
available.  To  fill  these  values  climatology  was  used.  The  sea  surface  temperature  (SST) 
and  air  temperature  were  variable  in  CHEMMAPTM  that  had  large  impacts  on  the 
outcome  of  the  Ammonia  and  Mustard  Gas  scenarios.  However,  climatologic  values  had 
to  be  used.  Eor  simplicity,  once  those  values  were  ingested  they  remained  the  same  for 
every  scenario.  Based  on  climatology,  a  SST  of  lO^C  was  utilized  (Walters  and  Sjoberg, 
1988).  In  the  northern  Persian  Gulf,  SST’s  are  almost  always  higher  than  this  value,  so  it 
allowed  for  less  evaporation  than  would  normally  occur. 

The  amount  used  in  the  scenarios  was  based  on  fields  obtained  from  HPAC. 
When  calculating  the  amount  of  Ammonia  to  be  released  in  the  HPAC,  several  options 
were  listed,  the  best  being  ‘industrial  transportation.  Of  those  options,  the  only  sea  option 
was  a  barge  that  contained  Liquid  Ammonia.  Since  Liquid  Ammonia  is  the  usual 
transportation  method  and  the  worst  case  scenario  made  the  best  scenario  to  study,  514 
tons  of  Liquid  Ammonia  entering  the  atmosphere  and  303  tons  being  pooled  was  utilized. 

Eor  Mustard  Gas  the  options  were  less  obvious.  When  articles  discussed  affects 
from  Mustard  Gas,  no  amounts  were  used.  Since  an  amount  was  needed,  the  same 
amount  used  for  the  Ammonia  scenarios  were  utilized.  This  would  give  an  extreme 
example  of  a  spill  scenario.  This  scenario  would  be  more  from  the  aspect  of  transporting 
the  chemical,  rather  than  an  actual  chemical  attack.  This  allowed  for  a  larger  amount  than 

any  chemical  attack  would  ever  likely  be  able  to  utilize.  Although  a  chemical  attack 
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could  be  in  the  form  of  an  airborne  attaek,  those  partieles  would  quiekly  land 
and  the  eollision  seenario  would  take  effeet. 

The  Mine  Warfare  (MIW)  Seenario  was  relatively  easy  to  design.  As  will  be 
discussed  later,  the  mines  were  ereated  to  have  a  1.0  m  diameter  and  density  of 
l.Qgjcrn  .  With  some  simple  caleulations,  it  was  determined  that  a  1.0  m  diameter  mine 
with  density  of  l.Oglcm^  would  equate  to  approximately  a  one  ton  mine.  Using  the  same 
amount  as  the  previous  seenarios  allowed  for  303  mines  to  be  laid  for  every  seenario. 


85 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


86 


IX.  MIW:  IMPACT  OF  DRIFTING  MINES  ON  OIL  PLATFORM 

OPERATIONS 


To  first  understand  the  results,  the  inputs  used  in  the  mine  drift  scenario  have  to 
be  examined.  The  scenario  started  with  a  generic  iron  filings  particle  in  CHEMMAP.  The 
particle  was  then  added  to  the  personal  database  for  modification.  To  simulate  an  actual 
mine  the  diameter  of  the  particle  was  enlarged  to  1.0  m.  Since  iron  is  much  denser  than 
water,  the  density  of  the  particle  had  to  be  modified  as  well.  To  allow  the  particle  to  float 
on  the  surface  and  simulate  a  possible  generic  floating  mine  the  density  was  changed 
tol.Og/cm^  .  The  scenario  also  included  inserting  numerous  mines  running  the  scenario 
and  then  changing  location  and  repeating  the  process.  Times  for  the  model  runs  were 
longer  for  MIW  scenarios  than  for  the  other  scenarios.  This  was  due  to  the  fact  that  once 
a  mine  is  placed  in  the  water  it  does  not  degrade  but  remains  until  it  is  either  hit,  washed 
ashore,  or  sinks.  As  such  scenarios  length  of  the  scenario  was  based  more  on  operational 
planning  than  on  anything  else. 

The  primary  scenario  that  was  examined  was  for  a  newly  laid  drifting  mine.  This 
scenario  could  be  examined  in  several  different  ways.  Some  of  these  could  be:  (1)  if  one 
mine  was  placed  in  the  water,  where  is  it  likely  to  travel;  (2)  if  a  large  amount  of  mines 
were  placed  in  the  water,  where  would  they  go  and  how  much  would  they  spread;  and 
lastly  (3)  given  that  a  mine  is  going  to  be  placed  in  the  water,  where  is  the  worst  case 
scenario  for  Iraq’s  oil  terminals,  or  is  there  a  place  that  the  oil  terminals  would  be 
affected. 

Drifting  mine  scenarios  were  run  using  different  locations  around  the  northern 
Persian  Gulf  (Figure  32).  As  can  be  seen  in  Figure  33,  when  the  mines  were  released 
from  locations  south  of  the  oil  terminals  the  mines  drifted  away  from  them  and  never 
affected  the  oil  terminal  operations  For  the  climatological  runs,  no  matter  where  the 
mines  were  placed,  they  followed  the  same  pattern  as  will  be  discussed  for  the  release 
near  the  oil  terminals.  This  was  primarily  due  to  the  nature  of  climatological  forcing. 
When  examining  scenarios  with  COAMPS  and  SWAPS  inputs,  the  worst  case  scenarios 
involved  mines  originating  near  the  oil  terminal.  This  is  because  mines  the  originated 
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south  of  the  oil  terminals  frequently  never  reached  the  oil  terminals  and  only  mines 
around  these  terminals  would  influence  operations  there.  Mines  located  away  from  these 
terminals  would  influence  traffic  coming  to  and  from  the  terminals  but  not  at  the 
terminals. 
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Figure  32.  Drop-off  points  utilized  during  the  drift  mine  scenario  model  runs.  White 

X’s  mark  mine  drop  locations,  while  the  two  stars  represent  oil  terminal  locations. 
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Figure  33.  Drifting  mine  scenario  with  release  from  a)  climatology  and  three  different 

release  locations  (b,  c,  and  d). 


A,  CLIMATOLOGICAL 

In  the  Persian  Gulf,  as  with  most  of  the  world,  there  is  no  climatological  data  for 
currents.  To  compound  matters,  there  were  no  current  meters  readily  available  to  verify 
the  physical  model  output.  Sporadic  current  observations  are  usually  accomplished 
during  set  research  studies.  Operationally,  this  presents  a  problem  when  in  the  planning 
phase.  What  would  the  planners’  tactics  be  to  account  for  currents  in  the  Persian  Gulf, 
when  none  are  available?  The  answer,  use  what  they  do  know,  tides.  In  all  likelihood 
either  a  tidal  model  or  a  cyclic  current  would  be  used  as  inputs  into  the  chemical  models 
for  planning  purposes.  With  this  concept  in  mind,  several  scenarios  were  used  around  the 
oil  terminals  and  northern  Persian  Gulf  to  simulate  planning  scenarios. 

The  first  scenario  utilized  a  cyclic  pattern  of  10  cm/s  to  simulate  a  tidal  pattern. 

This  pattern  did  what  was  expected.  The  mines  oscillated  around  the  point  where  they 
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were  dropped,  not  really  spreading  out  much.  When  the  climatological  wind  component, 
7  m/s,  was  added  the  mines  then  drifted  in  a  cluster  in  the  direction  of  the  climatological 
winds  (Figure  34,).  Of  note,  the  drift  would  follow  the  current  cycle.  There  were  periods 
of  faster  drift  with  the  wind,  followed  by  periods  of  slower,  sometimes  even  reversing 
direction,  when  the  current  opposed  the  climatological  winds. 


Figure  34.  Results  from  the  climatological  MIW  scenario,  which  included  cyclic 
currents  and  climatological  winds  over  a  5  day  period.  The  panels  are  for  a)  initialization, 

b)  day  2,  c)  day  4  and  d)  day  5. 

When  the  climatological  winds  were  used  without  any  currents,  the  mines  drifted 
in  the  direction  of  the  winds.  This  scenario  was  done  to  ensure  that  winds  were  indeed  a 
factor  on  the  simulated  mines.  Although  the  mines  drifted  with  the  climatological  winds, 
they  did  not  travel  at  the  same  speed  as  the  winds.  The  winds  input  were  7  m/s,  whereas 
the  mines  drifted  much  slower  than  this.  In  five  days  time  at  7  m/s  the  mines  would  have 
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drifted  over  600  km.  Sinee  the  Persian  Gulf  is  approximately  990  km  in  length,  if  the 
mines  were  traveling  at  the  same  speed  as  the  winds  they  would  have  traveled  well  over 
half  the  length  of  the  Persian  Gulf.  The  mines  however,  had  just  left  the  northern  region 
at  the  end  of  day  five. 

B,  SCENARIO  SPECIFIC 

Utilizing  modeled  data  was  the  next  step.  Several  periods  were  used  for  the  MIW 
seenario  and  subsequent  seenarios.  These  periods  were  speeifieally  ehosen  to  meet  one  of 
the  following  eriteria:  (1)  minimal  wind  and  weak  eurrent  event,  (2)  maximum  wind  and 
weak  current  event,  (3)  maximum  wind  and  maximum  current  event  and  (4)  looking  for 
some  period  that  a  weather  phenomenon  occurred.  The  weak  wind  and  low  current  period 
corresponded  to  20  through  25  February.  All  other  weak  wind  periods  either  did  not 
correspond  to  any  of  the  currents  criteria  or  was  too  close  to  bad  or  missing  data  to  be 
utilized.  During  this  event,  the  chemical  model  transports  the  mine  much  slower  than 
indicated  by  the  climatological  runs;  however  the  dispersion  is  greater  (Figure  35).  The 
most  likely  reason  for  greater  dispersion  is  the  high  resolution  of  currents  allowing  for 
more  dispersion.  With  a  uniform  field  in  the  climatological  scenario,  regardless  of  where 
the  mines  were  placed  they  felt  the  same  force.  Under  these  conditions,  little  dispersion 
would  be  expected.  In  February’s  scenario,  the  mines  felt  different  currents  depending  on 
where  they  were  located.  This  meant  if  the  starting  point  was  shifted,  mines  would 
experience  different  forcing  a  therefore  could  be  expected  to  disperse  more.  The  fact  the 
timeframe  examined  was  a  period  of  low  currents,  factored  in  to  the  amount  of 
dispersion.  As  the  other  scenarios  will  illustrate,  the  minimal  wind  and  weak  current 
scenario  had  the  least  dispersion  of  all  the  scenarios.  It  also  had  the  least  movement,  as 
would  be  expected. 
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c)  d) 

Figure  35.  Results  from  the  drifting  mine  scenario  during  a  weak  current  and  minimal 
winds  event  in  February  2006.  The  panels  depict  elapsed  time  of  a)  2  days,  b)  3  days,  c)  4 
days  and  d)  5  days.  Pink  dots  are  the  current  location  of  the  mines,  while  the  black  dots 

illustrate  track  history. 


There  were  two  events  that  met  the  low  currents  and  maximum  winds  criteria. 
The  first  was  from  7  through  12  April.  In  the  April  scenario  the  winds  frequently  met  the 
high  wind  criteria,  but  would  also  repeatedly  decrease  below  the  threshold.  The  other 
event  occurred  between  3  and  10  July.  The  July  event  was  significant  because  the  winds 
were  the  highest  seen  during  the  entire  period,  while  the  currents  continued  to  stay  below 
the  25  cm/s  criteria.  What  is  observed  from  running  the  chemical  model  during  these 
periods  was  significantly  different  from  the  climatological  outputs  from  the  previous 
section.  Both  periods  depicted  the  mines  dispersing  further  than  climatology  would 
suggest.  When  examining  the  results  from  the  April  scenario.  Figure  36,  the  mines  can  be 
seen  to  disperse  more  than  the  climatologic  scenario;  however  they  did  not  travel  as  far. 
This  is  likely  due  to  the  winds  in  the  scenario  being  less  than  what  climatology  suggests. 
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The  winds  in  the  July  event  are  eloser  to  the  elimatologieal  means  and  traveled  much 
further.  The  July  scenario  also  saw  a  large  amount  of  dispersion  (Figure  37). 
Concentrations  of  mines  spread  in  all  direction,  yet  as  the  chemical  model  run  progressed 
an  East-West  line  reformed  (Figure  37).  The  larger  wind  speeds  likely  contributed  to  the 
increased  movement.  However  the  eventual  East-West  spreading  was  an  unexpected 
result.  All  other  model  runs  tended  to  leave  more  of  a  circular  or  elliptical  pattern.  A 
possible  reason  for  the  outcome  was  the  length  of  the  run.  All  other  model  runs  for  the 
drifting  mine  scenario  were  for  five  days,  the  July  event  was  run  for  ten  days.  The  output 
at  the  five  day  point.  Figure  37  was  similar  to  the  other  events.  A  further  study  utilizing 
longer  runs  may  be  beneficial  for  operational  planning. 


c)  d) 

Figure  36.  Results  from  the  drifting  mine  scenario  during  a  weak  current  and 
maximum  winds  event  in  April  2006.  The  panels  depict  elapsed  time  of  a)  1  day,  b)  3 

days,  c)  4  days  and  d)  5  days. 
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The  criterion  for  high  currents  combined  with  high  winds  was  seen  once  from  9 
until  16  June.  Although  the  high  current  criterion  was  met,  the  currents  oscillated 
direction  at  high  speeds  and  did  not  maintain  one  direction  (Figure  38).  This  led  to  a 
shorter  distance  traveled  with  the  same  dispersion  seen  in  the  July  scenario  (Figure  37). 
This  includes  forming  a  line  similar  to  the  July  event.  The  main  difference  is  that  the  line 


would  form,  break  apart  and  then  reform. 


c)  d) 

Figure  37.  Results  from  the  drifting  mine  scenario  during  a  weak  current  and 
maximum  winds  event  in  July  2006.  The  panels  depict  elapsed  time  of  a)  1  day,  b)  3 

days,  c)  4  days  and  d)  5  days. 


The  primary  conclusion  from  the  drifting  mine  scenario  is  that  regardless  of 
where  mines  are  placed  in  the  northern  Persian  Gulf  they  tend  to  move  away  from  the  oil 
terminals.  This  fact  means  that  while  ships  transiting  to  and  from  the  terminals  are 
affected,  those  around  the  terminals  would  not  be  affected.  This  result  implies  minimal 
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impact  to  the  actual  operation  of  the  oil  terminals;  however,  if  ships  eould  not  reach  the 
terminals  operations  would  be  limited. 


e)  d) 


Figure  38.  Results  from  the  drifting  mine  seenario  during  a  strong  oseillating  eurrent 
and  maximum  winds  event  in  June  2006.  The  panels  depict  elapsed  time  of  a)  2  days,  b) 

3  days,  c)  4  days  and  d)  5  days. 


95 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


96 


X.  IMPACT  OF  A  CHEMICAL  SPILL  TO  OPERATIONS: 

AMMONIA 


The  Ammonia  scenario  included  both  an  atmospheric  portion  and  oceanographic 
portion.  The  Atmospheric  portion  involved  running  HP  AC,  while  the  oceanographic 
portion  involved  CHEMMAP™.  Although  CHEMMAP™  contains  an  atmospheric 
model,  the  primary  feature  examined  was  oceanographic. 

A,  ATMOSPHERIC  PORTION  OF  SCENARIO 

The  atmospheric  portion  of  the  Ammonia  scenario  drove  the  amounts  utilized  for 
the  oceanographic  portion.  Amounts  were  determined  through  release  type.  Since  the 
goal  was  to  look  at  a  worst  case,  yet  realistic,  scenario,  it  was  determined  that  a 
transportation  scenario  was  the  best  choice.  This  option  allowed  for  total  damage  of  a 
Type  3  Barge  carrying  Ammonia.  The  scenario  also  automatically  calculated  the  amount 
of  the  spill  that  entered  the  water  and  the  water.  The  total  viable  airborne  Ammonia  was 
calculated  to  be  514  tons.  Of  this  amount,  258  tons  were  in  vapor  form,  while  the 
remaining  257  tons  were  an  aerosol.  HPAC  also  calculated  that  303  tons  of  Ammonia 
would  be  pooled  in  the  water.  The  only  changes  that  could  be  made  to  these  values  were 
to  change  the  damage  amount  or  the  transportation  type.  Since  the  scenario  was  looking 
at  worst  case,  total  damage  was  most  appropriate.  With  the  intent  of  the  scenario  to  affect 
Iraq’s  oil  terminals,  the  only  viable  option  was  a  waterborne  vehicle.  Within  HPAC  this 
left  only  a  Type  3  Barge.  Since  the  amount  within  HPAC  was  extremely  rigid,  and  based 
on  genuine  values,  the  scenario  was  utilized  for  the  oceanographic  portion  of  the 
scenario. 

1,  Climatological 

One  of  the  features  imbedded  in  HPAC  is  the  ability  to  choose  climatology  for  a 

given  period  rather  than  real-time  or  archived  data.  This  climatological  data  allows  the 

user  to  choose  the  month  and  that  field  is  used  for  the  model  run.  Some  problems  with 

this  include:  (1)  what  month  to  choose  if  the  data  covers  the  end  of  one  month  and  the 

beginning  of  the  next,  and  (2)  trying  to  determine  the  values  for  the  climatological  wind 

field  that  is  being  used.  The  first  issue  was  avoided  since  the  times  chosen  never 

overlapped  months.  This  was  primarily  done  to  allow  currents  data  to  be  utilized  within 
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CHEMAP^"^,  however  the  desired  criteria  also  occurred  in  such  a  manner  as  to  avoid  this 
issue  altogether.  Trying  to  extract  the  climatological  winds  proved  futile.  Although  from 
the  outputs  an  important  conclusion  could  be  made.  The  climatological  winds  were 
constant  in  both  speed  and  direction. 

As  with  the  climatological  MIW  scenario  within  CHEMMAP™,  all 
climatological  times  resulted  in  the  same  outcomes.  Therefore,  only  one  scenario  will  be 
illustrated.  This  scenario,  top  half  of  Eigure  39,  produced  a  cone-like  dosage  duration 
plume.  This  plume  extended  from  near  the  oil  terminals  along  the  major  axis  of  the 
Persian  Gulf.  Although  this  plume  looks  significant,  the  maximum  time  for  duration  was 
ten  minutes.  With  this  in  mind,  climatology  indicated  that  a  chemical  spill  would  rapidly 
leave  the  area  and  transit  south.  Other  methods  of  looking  at  the  same  scenario 
demonstrated  chemical  would  also  spread  out  quite  rapidly.  Within  a  days  time  the 
concentration  had  migrated  away  from  the  sight  of  interest  and  had  spread  so  that  the 
concentrations  were  not  harmful  to  humans  (lower  half  of  Eigure  49). 
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c)  d) 

Figure  39.  The  upper  portion  depicts  a)  50  minutes  and  b)  15  hours  elapsed  of  514 
tons  of  Ammonia  for  the  integrate  concentration.  The  lower  portion  indicates  c)  50 
minutes  and  d)  18  hours  elapsed  for  the  duration  of  Ammonia  in  the  atmosphere. 


2,  Specific  Scenario 

Utilizing  modeled  data  was  the  next  step.  The  same  criteria  utilized  in  the  MIW 
scenario  were  initially  utilized  to  determine  model  times  periods  for  the  atmospheric 
portion  of  the  Ammonia  scenario.  In  the  end  the  atmospheric  criteria  were  changed 
slightly.  The  criteria  now  became:  (1)  relatively  constant  onshore  flow,  (2)  relatively 
constant  offshore  flow,  and  (3)  include  an  event  such  as  frontal  passage  in  the  analysis. 
Another  difference  was  the  timescale.  Since  the  climatological  runs  demonstrated  that 
Ammonia  would  leave  the  zone  of  interest  much  faster  than  in  the  MIW  scenario,  a 
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shorter  model  run  was  utilized.  Therefore,  the  model  was  run  over  a  24  hour  period.  This 
shorter  period  allowed  for  more  model  runs.  Since  there  were  not  large  differences  in  the 
results,  only  the  worst  case  scenario  will  be  shown  for  each  criterion. 

Running  the  model  in  late  March  met  the  criteria  of  relatively  constant  onshore 
flow.  The  analysis  of  this  criterion  also  included  several  times  where  frontal  patterns 
were  observed.  These  runs  included  a  wind  shift  from  an  onshore  flow  to  an  alongshore 
flow.  By  the  time  the  flow  had  shifted  the  scenario  was  complete,  however  the  criteria  for 
offshore  flow  would  soon  be  met.  A  problem  encountered  when  doing  this  was  onshore 
flow  rapidly  approached,  and  left,  the  archived  domain.  Once  this  happened,  the  fate  of 
the  chemical  could  only  be  known  to  have  left  the  modeled  region.  From  Figure  40,  it 
was  determined  that  after  15  hours  the  majority  of  the  chemical  had  traveled  outside  the 
modeled  domain.  It  was  also  determined  that  after  the  first  four  hours  levels  were  below 
threshold  limits  from  ammonia,  with  possible  injury  no  longer  being  a  factor  after  two 
hours. 
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a)  b) 

Figure  40.  Trajectory  of  5 14  tons  of  Ammonia  at  elapsed  times  of  a)  6  hours  and  b)  15 
hours  during  an  onshore  winds  scenario  in  March  2006. 

The  other  scenario  involved  offshore  flow.  This  scenario  occurred  more  often  and 
was  more  inline  with  climatology.  As  seen  in  Figure  41,  the  offshore  flow  rapidly 
dispersed  the  chemical  throughout  the  region.  This  also  meant  that  the  toxic  levels 
dropped  off  almost  immediately.  From  the  results  of  these  two  scenarios  it  can  be 
concluded  that  even  a  major  Ammonia  spill  will  not  have  very  long  lasting  effects  for 
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humans.  Other  studies,  sueh  as  Claereboudt,  et  al.  (2001),  have  suggested  that  higher 
Ammonia  levels  will  drastieally  affect  aquatic  life  in  the  region.  However,  this  study 
focused  on  operational  impacts,  which  lasted  much  less  than  a  day. 


c)  d) 

Figure  41 .  Trajectory  of  5 14  tons  of  Ammonia  at  elapsed  times  of  a)  6  hours,  b)  12 

hours,  c)  18  hours  and  d)  24  hours  during  an  offshore  winds  scenario  in  May  2006. 


B.  OCEANOGRAPHIC  PORTION  OF  SCENARIO 

1,  Climatological  Input 

The  climatology  utilized  for  the  currents  is  the  same  as  in  the  MIW  scenario. 
Since  no  real  climatology  exists,  the  climatology  utilized  was  in  reference  to  what 
planners  may  consider  during  the  planning  phase.  Not  surprisingly,  the  results  tended  to 
mirror  that  of  the  MIW  scenario.  One  of  the  only  differences  was  in  the  time  scale.  Since 
mines  do  not  degrade  or  evaporate  a  longer  timescale  was  used  for  them,  whereas  for  the 
Ammonia  scenario  degradation  and  evaporation  had  to  be  taken  into  account  when 
running  the  scenario.  After  running  several  climatologically  forced  scenarios,  it  was 
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determined  that  the  Ammonia  deeayed  on  the  hour  time  frame,  not  the  day  time  scale 
(Figure  42).  Running  the  models  for  48  hours  produced  the  best  results.  After  48  hours 
the  concentration  was  barely  detectable.  The  scenario  also  produced  a  narrow  plume  that 
decreased  in  intensity  over  time.  The  climatology  scenario  moved  the  Ammonia  farther 
away  and  in  a  straight  line.  This  was  due  to  the  uniform  forcing  and  lack  of  variability 
acknowledged  previously  in  the  mine  drift  scenario. 


a)  b) 

Figure  42.  Two  panel  view  depicting  a)  24  hours  and  b)  72  hours  elapsed  time  of 

Ammonia  with  climatological  forcing. 


2,  Specific  Scenario 

Utilizing  modeled  data  was  the  next  step.  While  the  climatological  scenario  was 
run  for  each  month,  the  time  periods  for  the  archived  data  utilized  the  same  philosophy  as 
in  the  MIW  scenario.  The  only  difference  was  the  timescale.  As  discussed  earlier  in  this 
chapter,  the  influences  of  degradation  and  evaporation  had  to  be  taken  into  account  when 
running  the  scenario.  Furthermore,  the  time  scale  should  match  the  runs  utilized  in  the 
climatology  portion.  Therefore,  48  hours  was  used  as  the  primary  length.  As  can  be  seen 
in  Figure  43,  the  majority  of  the  liquefied  Ammonia  evaporated  within  3  hours.  Whatever 
did  not  evaporate  quickly  spread  throughout  the  water  column.  The  fact  that  the  water 
column  was  around  30  m  thick  allowed  this  concentration  to  continually  evaporate  over 
the  next  12  hours.  There  was  also  a  small  percent  of  Ammonia  that  decayed.  The 
concentrations  left  after  12  hours  remained  near  10  percent  of  the  total  concentration 
spilled.  This  amount  is  the  concentration  depicted  in  Figure  44.  The  figure  clearly  shows 
low  concentrations  in  the  water  and  very  little  flow.  The  limited  movement  of  Ammonia 
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was  obtained  due  to  eurrents  being  primary  influence  on  Ammonia.  Since  the  Ammonia 
quickly  entered  the  water  column,  wind  had  little  effect.  Although  the  amount  of  the  spill 
was  designed  to  be  worst  a  worst  case  scenario,  the  impact  on  the  oil  terminals  were 
minimal.  The  fact  is  that  although  Ammonia  is  abundant  in  the  region,  alone  it  does  not 
pose  a  long  term  threat  to  Iraq’s  oil  terminal  operations  with  in  the  Persian  Gulf 

Mass  Balance  for  Ammonia-liquefied  gas 


Jul  2006  Date/Time 

Figure  43.  12  hour  mass  balance  graph  of  Liquefied  Ammonia  gas  within 

CHEMMAP™. 
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Figure  44.  Four  panel  view  depieting  a)  initialization,  b)  72  hours  elapsed  during  April 
2006,  e)  72  hours  elapsed  time  during  June  2006and  d)  72  hours  elapsed  time  during  July 
2006  for  the  Ammonia  spill  scenario.  The  colorbar  units  are  in  mg/  m^ .  The  labels  are  as 
follows;  Tan  =  0.01-0.1,  Yellow  =  0. 1-0.5,  Light  Green  =  0. 5-1.0,  Green  =  1-5,  Light 
Blue  =  5-10,  Blue  =  10-50,  Light  Purple  =  50-100,  Purple  =  100-500,  and  Orange  >  500. 
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XI.  IMPACT  OF  A  CHEMICAL  SPILL  TO  OPERATIONS: 

MUSTARD  GAS 

The  biggest  influenee  in  the  Mustard  Gas  scenario  was  SST.  With  water 
temperatures  being  unavailable  from  the  archived  fields,  climatological  fields  had  to  be 
used  for  all  scenarios.  This  meant  that  no  low  values,  which  would  have  slowed  the 
degradation,  were  available.  It  also  meant  that  to  get  a  worst  case  scenario,  values  lower 
than  climatology  would  be  used.  In  CHEMAP™,  values  for  water  temperature  would  be 
set  at  20°C.  This  value  is  lower  than  climatology  suggests  and  therefore  would  be 
considered  a  worst  case  scenario  for  evaporation  and  degradation. 

The  high  SST  turned  out  to  be  the  biggest  factor  in  all  Mustard  Gas  scenarios. 
Regardless  of  other  inputs,  the  high  SST  combined  with  high  air  temperatures  meant 
model  runs  were  limited  to  a  maximum  of  12  hours  (Figure  45).  This  figure  was  similar 
for  every  model  run  regardless  of  meteorological  and  oceanographic  forcing  or  the  time 
when  it  was  run.  To  prove  that  temperature  was  indeed  a  major  influence,  several  model 
runs  were  conducted  using  30°C  and  10°C.  Although  these  values  are  not  realistic  in  the 
northern  Persian  Gulf,  they  were  valuable  for  proving  the  influence  of  SST  on  Mustard 
Gas  duration.  As  Figures  46  illustrates,  the  colder  temperature  decreased  the  degradation 
rate  and  increased  the  time  Mustard  Gas  remained  in  the  area,  whereas  Figure  47 
demonstrated  the  opposite  effect.  From  Spectrum  Faboratories  (2007)  and  Floyd  (2007) 
most  of  the  effects  from  Mustard  Gas  were  from  evaporation  during  World  War  II.  The 
problem  was  that  the  chemical  attacks  would  occur  at  night  and  then  the  chemicals  would 
evaporate  during  the  day.  This  is  less  likely  to  happen  in  the  northern  Persian  Gulf  due  to 
the  higher  SST’s  and  air  temperature.  With  this  in  mind  and  the  fact  the  scenarios  run 
were  for  a  major  collision,  it  can  be  concluded  that  the  initial  impact  to  a  Mustard  Gas 
attack  would  not  effect  oil  terminal  operations  longer  than  two  hours.  The  additional 
delay  would  be  due  to  clean-up,  removal  of  injured  personnel  and  bringing  new  personnel 
onboard.  The  additional  time  for  clean-up  would  be  minimal  due  to  the  rapid  decay  seen 
in  Figure  45.  The  entire  time  from  initial  attack  to  a  full  operational  oil  terminal  would  be 
measured  in  hours.  Therefore,  the  result  of  this  study  is  that  utilizing  any  other  data 


105 


besides  temperature  does  not  have  any  perceivable  value  added.  This  is  due  to  the  short 
duration,  which  leads  to  minimal  impact  from  other  fields  in  the  environment. 
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Figure  45.  6  hour  mass  balance  graph  of  Mustard  Gas  at  20“C  within  CHEMMAP™. 
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Figure  46.  6  hour  mass  balance  graph  of  Mustard  Gas  at  30“C  within  CHEMMAP™. 
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XII.  SUMMARY 


Starting  with  six  months  of  surface  winds  and  currents  data  obtained  from 
CO  AMPS  and  SWAPS,  an  analysis  was  completed  to  determine  what  the  value  added 
would  be  to  a  drifting  mine  scenario  and  two  different  chemical  spill  scenarios  in  the 
northern  Persian  Gulf.  This  region  was  determined  to  be  a  low-energy  region.  The 
region’s  primary  wind  influences  were  determined  to  be  the  Shamal  and  transiting  low 
pressure  systems.  The  primary  influence  of  the  currents  corresponded  to  wind-driven 
currents  and  tides. 

The  drifting  mine  scenario  included  several  hundred  simulated  mines  placed  in 
various  locations  throughout  the  northern  Persian  Gulf  region.  The  worst  case  scenario 
was  determined  to  be  when  mines  were  placed  north  of  the  oil  terminals  or  near  them. 
This  scenario  gave  the  best  chance  of  mines  to  affect  the  oil  terminals.  If  the  mines  were 
placed  to  the  south,  the  current  would  normally  carry  them  away.  The  mine  scenario  had 
significant  value  added  using  any  real-time  data,  as  opposed  to  trying  to  use  climatology. 
It  was  also  determined  that  climatological  data  is  not  available  for  currents  and  that, 
during  a  planning  phase,  either  cyclic  or  tidal  currents  would  be  used.  These  currents  did 
not  disperse  the  mines  as  much  as  scenarios  involving  the  archived  data.  In  the  end,  even 
if  a  weak  event  was  utilized,  the  dispersion  was  greater  than  climatological  methods,  and 
therefore  more  realistic. 

The  Ammonia  scenario  contained  an  entirely  atmospheric  portion  as  well  as  an 
oceanographic  portion.  This  scenario  utilized  a  collision  scenario  to  provide  the  worst 
case  for  an  Ammonia  spill.  The  results  clearly  demonstrated  that  using  real-time,  or  for 
this  study  archived,  model  data  would  drastically  change  the  results  from  using 
climatological  data.  In  the  atmospheric  portion,  two  events  demonstrated  a  large 
discrepancy  between  climatology  and  modeled  data.  The  event  most  like  climatology,  the 
offshore  event,  produced  a  greater  dispersion  of  the  Ammonia  and  allowed  it  to  travel 
much  farther  as  well.  The  other  event  had  the  flow  in  the  opposite  direction  to 
climatology.  This  onshore  event  occurred  primarily  when  transiting  lows  came  through 

the  area.  Although  this  event  produced  results  drastically  different  from  climatology,  the 

109 


results  were  limited  by  the  Ammonia  rapidly  leaving  the  modeled  area.  Operational 
planners  would  have  to  create  two  scenarios.  These  two  scenarios  would  cover  the  worst 
case  events  for  the  period  analyzed.  These  events  would  be:  (1)  an  offshore  event  that 
would  quickly  spread  the  Ammonia,  and  (2)  an  onshore  event. 

The  other  portion  of  the  Ammonia  scenario  was  conducted  with  in 
CHEMMAP^'^.  In  this  scenario  the  same  spill  amount  was  used  but  the  oceanographic 
portion,  the  pooled  303  tons,  was  analyzed.  The  majority  of  the  liquefied  Ammonia 
evaporated  within  three  hours.  Therefore  it  was  determined  that  Ammonia  did  not  stay  in 
the  water  long  enough  to  have  much  of  an  operational  impact,  and  that  aquatic  life  would 
be  the  primary  subject  effected. 

The  last  scenario  analyzed  was  a  chemical  spill  of  Mustard  Gas.  This  event  used 
the  same  amount  as  in  the  Ammonia  scenario.  Although  this  amount  is  much  larger  than 
would  be  used  in  a  chemical  attack  scenario,  it  could  occur  if  there  was  a  collision  when 
chemical  weapons  were  being  moved.  Although  303  tons  of  Mustard  Gas  were  released, 
within  one  to  two  hours  all  the  chemical  had  evaporated  or  decayed.  This  was  due  to  the 
high  SSTs  and  air  temperatures  found  in  the  northern  Persian  Gulf.  These  high 
temperatures  were  determined  to  be  the  limiting  factors  for  Mustard  Gas.  Since  the  data 
collected  only  included  wind  fields,  no  temperature  data  was  available  for  further 
analysis.  The  temperature  fields  that  were  used  during  this  study  were  kept  below  the 
climatological  means  to  provide  worst  case  scenarios.  Even  by  using  these  lower  values, 
there  was  little  impact,  beyond  the  initial  strike/collision  and  those  immediately  effected, 
to  oil  terminal  operations.  The  fact  the  operations  would  most  likely  be  able  to  resume 
within  a  couple  of  hours  meant  that  no  value  was  added  by  using  any  wind  values,  either 
climatology  or  modeled. 

As  regards  recommendations  for  future  research,  it  should  be  mentioned  that  the 
mine  drift  scenarios  utilized  CHEMMAP^'^  to  generate  a  generic  mine.  This  mine  was 
1.0  meters  in  diameter  and  had  a  density  of  l.Qgjcrn  .  Since  this  size  and  density  was 
not  based  on  any  real  data,  further  study  utilizing  CHEMMAP^"^  with  more  realistic 
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floating  mine  density  and  sizes  could  prove  beneficial.  If  ASA  could  add  a  drifting-mine 
database  to  CHEMMAP^'^,  the  program  could  be  useful  for  the  MIW  operations. 

If  archived  temperature  fields  could  be  retrieved  for  the  given  timeframe,  a 
detailed  analysis  of  those  fields  could  give  a  better  picture  of  temperature  ranges  in  the 
northern  Persian  Gulf.  These  ranges  would  assist  in  determining  the  exact  duration  of 
Mustard  Gas  and  then  if  there  are  any  impacts  on  operations. 
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APPENDIX 


This  appendix  contains  all  relevant  figures  not  referred  to  in  the  text. 
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Figure  48.  EOF-1  for  winds  of  the  entire  Persian  Gulf 
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SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  Model) 
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Figure  49.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  wind  dataset  in  the 

Persian  Gulf 
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Figure  50.  EOF-2  for  winds  of  the  entire  Persian  Gulf 
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Figure  5 1 .  EOF-2  amplitude  (top)  and  direction  (bottom)  for  entire  wind  dataset  in  the 

Persian  Gulf 
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Figure  52.  EOF-3  for  winds  of  the  entire  Persian  Gulf 


SWAPS  Magnitude/Theta  Time  Series:  (Magnitude  Mode3) 


Figure  53.  EOF-3  amplitude  (top)  and  direetion  (bottom)  for  entire  wind  dataset  in  the 

Persian  Gulf 
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Figure  54.  EOF-1  for  currents  of  the  entire  Persian  Gulf.  Vectors  are  every  fourth 
vector  of  the  ones  calculated,  which  equates  to  every  eighth  point. 
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Figure  55.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  current  dataset  in 

the  Persian  Gulf. 
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Figure  56.  EOF-2  for  currents  of  the  entire  Persian  Gulf.  Vectors  are  every  fourth 
vector  of  the  ones  calculated,  which  equates  to  every  eighth  point. 

SWAPS  Magnitude/Theta  Time  Series  (Magnitude  Mode  2) 


o'  '  '  '  . .  " 

1  Feb  1  Mar  29  Mar  26  Apr  24  May  21  Jun  19Jul 


(Theta  Mode  2) 


1  Feb  1  Mar  29  Mar  26  Apr  24  May  21  Jun  19Jul 

Figure  57.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  current  dataset  in 

the  Persian  Gulf. 
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Figure  58.  EOF-3  for  currents  of  the  entire  Persian  Gulf.  Vectors  are  every  fourth 
vector  of  the  ones  calculated,  which  equates  to  every  eighth  point. 
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Figure  59.  EOF-1  amplitude  (top)  and  direction  (bottom)  for  entire  current  dataset  in 

the  Persian  Gulf. 
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